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In this paper, we perform a detailed analysis of the Exact Holographic Mapping first introduced in 
arXiv: 1309.6282, which was proposed as an explicit example of holographic duality between quantum 
many-body systems and gravitational theories. We obtain analytic results for free fermion systems 
that not only confirm previous numerical results, but also elucidate the exact relationships between 
the various physical properties of the bulk and boundary systems. These analytic results allow us 
to study the asymptotic properties that are difficult to probe numerically, such as the near-horizon 
regime of the black hole geometry. We shall also explore a few interesting but hitherto unexplored 
bulk geometries, such as that corresponding to a boundary critical fermion with nontrivial dynamical 
critical exponent. Our analytic framework also allows us to study the holographic mapping of some 
of these boundary theories in dimensions 2+1 or higher. 
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I. INTRODUCTION 

In the recent years, holographic duality, also known 
as the Anti-de-Sitter space/Conformal Field Theory 
(AdS/CFT) correspondence[1-3], has attracted tremen¬ 
dous research interest in both high energy and condensed 
matter physics. This correspondence is defined as a du¬ 
ality between a D -\- 1-dimensional field theory on a fixed 
background geometry and a D + 2-dimensional quan¬ 
tum gravity theory. The best understood example of 
holographic duality is the correspondence between 4- 
dimensional super-Yang-Mills theory and 5-dimensional 
supergravity. There, the large-7V limit of the super- 
Yang-Mills theory corresponds to the classical limit of the 
dual gravity theory, which provides a helpful description 
of strongly coupled gauge theories. What makes holo¬ 
graphic duality particularly interesting is its generality. 
When the boundary theory is not a conformal field the¬ 
ory, a dual theory with a different space-time geometry 
may still be well-defined. [4] Physically, holographic du¬ 
ality can be understood as a generalization of the renor¬ 
malization group (RG) flow of the boundary theory [5-7], 
where bulk gravitational dynamics generalize the RG flow 
equations and the emergent dimension perpendicular to 
the boundary has the physical interpretation of energy 
scale[8, 9]. Indeed, holographic duality has been applied 
to condensed matter physics as a new tool to characterize 
strongly correlated systems [10-13]. 

More recently, holographic duality has been proposed 
to be related to another approach developed in con¬ 
densed matter physics, namely tensor networks[14-22]. 
In its most general form, tensor networks refer to a 
description of many-body wavefunctions and operators 
(i.e. linear maps) by contracting tensors defined on ver¬ 
tices of a graph. [23-28] More specifically, the tensor net¬ 
work state proposed to be related to holographic dual¬ 
ity is the multiscale entanglement renormalization ansatz 
(MERA)[26, 27], which is defined on a graph with hyper¬ 
bolic structure, with external indices (corresponding to 
the physical degrees of freedom) at the boundary and 
internal indices contracted in the bulk. An important 


feature of states described by tensor networks is that 
the entanglement entropy of a given region is bounded 
by the number of links between the region and its com¬ 
plement. This property motivated its relation to holo¬ 
graphic duality [14], where the entanglement entropy of a 
given region is determined by the area of the minimal sur¬ 
face bounding it, in accordance to the Ryu-Takayanagi 
formula [29]. 

There are many open questions in the proposed tensor 
network interpretation of the holographic duality. One 
important question is how to describe space-time geome¬ 
try rather than spatial geometry. Another (related) ques¬ 
tion is how to understand excitations (quantum fields) 
living in the bulk. Motivated by these questions, one 
of us [30] proposed a tensor network which defines not 
a many-body state but a unitary mapping between the 
boundary and bulk systems, known as the exact holo¬ 
graphic mapping EHM. The EHM is a tensor network 
very similar to MERA, except that it is a one-to-one uni¬ 
tary mapping between boundary and bulk degrees of free¬ 
dom. Each boundary state | ip) is mapped to a bulk state 
|^) = M |^), and each boundary operator O is mapped to 
a bulk operator O = MOM -1 . Physically, the EHM is 
a “lossless” version of real space renormalization group. 
Denoting a site in the bulk as x, a local operator at that 
site O x is dual to a generically nonlocal operator on the 
boundary O x = M _1 O x M. Different bulk sites x corre¬ 
spond to operators O x on the boundary with different en¬ 
ergy scales and different center-of-mass locations of their 
support. Once a mapping M is chosen, bulk correla¬ 
tion functions can in principle be calculated. Motivated 
by the general principle of relativity, the bulk geometry 
was proposed to be determined by the bulk correlation 
functions. More specifically, the distance between two 
points was proposed to depend logarithmically on the 
connected two point correlation functions. Compared to 
previous tensor network proposals, the EHM is different 
in two aspects: i) The bulk geometry is not determined 
by the structure of the tensor network but by the corre¬ 
lation structure of the bulk state; ii) The bulk geometry 
can be studied in both the spatial and temporal direction 
by studying the bulk correlation functions. In Ref. 30, 
an explicit choice of the mapping M for (1 + l)-dim lat¬ 
tice fermions was proposed, and the consequent dual ge¬ 
ometries corresponding to different boundary states were 
studied. They included the ground state of massless and 
massive fermions, the nonzero temperature thermal en¬ 
semble of massless fermions, and a thermal double state 
which is a purification of the thermal ensemble. Dy¬ 
namics after a quantum quench was also studied in the 
thermal double system, motivated by a comparison with 
geometrical properties of a two-sided black hole space- 
time[18]. 

The results in Ref. 30 for the abovementioned free 
fermion systems were obtained numerically. This limits 
the extent of analysis, due especially to the exponential 
growth of boundary system size. To have a well-defined 
bulk geometry with N layers in the emergent direction 
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of the bulk perpendicular to the boundary, the boundary 
system has to have 2 N sites. In this paper, we shall ob¬ 
tain analytic results on the free fermion EHM, which will 
enable us to rigorously determine asymptotic properties 
of the dual geometry, and also to discuss more general 
boundary systems. For instance, the existence of a black 
hole horizon in the geometry dual to a nonzero temper¬ 
ature state at the boundary can be studied more explic¬ 
itly from the asympotic infrared behavior of correlation 
functions in both spatial and temporal directions. In ad¬ 
dition to reproducing the results of Ref. 30 analytically, 
we shall also explore a few other interesting emergent 
bulk geometries, such as that corresponding to a critical 
fermion with nontrivial dynamic critical exponent. Our 
analytic framework also allows us to generalize the EHM 
to boundary theories with dimension 2 + 1 or higher [31], 
in which case the analytic approach is more essential due 
to the increasing difficulty of numerical calculations [32]. 
An added advantage of an analytic approach is that it al¬ 
lows one to identify properties of the bulk geometry that 
are insensitive to details of the choice of the mapping 
which thus reflects intrinsic properties of the boundary 
state. 

This paper is structured as follows. In Section II, we 
first review the EHM construction by describing its gen¬ 
eral principles and the definition of bulk geometry. These 
ideas will be elaborated in Section III for free lattice 
fermions, where an explicit Haar wavelet representation 
of the EHM will be presented. In Section IV, we provide 
detailed descriptions of the asymptotic correlator behav¬ 
ior and corresponding bulk geometries for the prototyp¬ 
ical 1 + 1-dim Dirac model at various combinations of 
zero and nonzero temperature and mass. These develop¬ 
ments will be further extended to higher dimensions and 
generic energy dispersions in Section V, where we discuss 
the emergence of interesting geometries like anisotropic 
black hole horizons with nontrivial topology. 

II. REVIEW OF THE EXACT HOLOGRAPHIC 
MAPPING 

In this section, we shall review the motivation and con¬ 
struction of the EHM proposed in Ref. 30 in a formalism 
that will be helpful for the later part of this paper. We 
will also include some new insights that are not discussed 
in the original proposal. The EHM approach is defined 
by the following two principles: 

1. The bulk theory and boundary theory are defined 
in the same Hilbert space. The bulk local operators 
are determined by a unitary mapping acting on the 
boundary local operators. 

2. The bulk geometry is determined by physical cor¬ 
relation functions. More specifically, the distance 
between two space-time points x, y in the bulk is 
determined by the connected correlation functions 
between the two points. 


Although the abovementioned unitary transformation 
can be very generic in principle, the types of transforma¬ 
tions that are relevant for holographic duality are those 
which are physically analogous to the renormalization 
group[33, 34]. The bulk operators at different locations 
should represent boundary degrees of freedom with dif¬ 
ferent energy scales. The key difference from the con¬ 
ventional RG approach is that the high energy degrees 
of freedom are spatially separated from low energy ones, 
instead of being integrated out. This enables us to con¬ 
cretely answer many new questions, such as how the high 
and low energy degrees of freedom (DOFs) are entan¬ 
gled/correlated. In the following, we will elaborate on 
the two abovementioned principles in the context of free 
fermion systems, and discuss the transformation of free 
fermion Hamiltonians under EHM. 


A. General construction of EHM 

The Exact Holographic Mapping is a unitary trans¬ 
formation defined by a tensor network or, equivalently, 
a quantum circuit consisting of local unitary operators. 
As proposed in Ref. 30, a simple construction of the 
EHM is given by a tree-shaped tensor network depicted 
in Fig. 1, where bulk (red) sites at the same level be¬ 
long to the same ’layer’. To construct it, we first take 
a D + 1-dimensional boundary system to be the zeroth 
bulk layer with L D = 2 nd sites. To construct the first 
bulk layer, one performs a unitary transform U on ev¬ 
ery set of 2 d adjacent sites such that the UV and IR 
(high and low momentum, assuming a monotonic energy 
dispersion) degrees of freedom are separated out. For 
D = 1, this can be written as 

Ul2^ lV>2> (1) 

a, (3 

where U \2 only acts on states on sites 1 and 2 re¬ 

spectively, and | a) and \/3) capture the higher and lower 
momentum (shorter and longer scale) degrees of freedom 
respectively. The construction of these \a) and \/3) states 
will be shown in detail in the next section. The full trans¬ 
formation on the zeroth layer is given by 

U = U 12 G C /34 G ... G U 2 N_i^ N , ( 2 ) 

which is a unitary transform on the Hilbert space of the 
whole layer. For D > 1 dimensions, U will be given by 
the direct product of D copies of the expression in Eq. 
2 . 

We construct the first bulk layer from the component 
|/?) = \/3)uv in Eq. 1 , which has the UV half of the 
degrees of freedom in the original layer. The other lower 
energy half |a), which we shall call the auxiliary sites 
in deference to Ref. 30, are fed into another copy of U 
with half the number of sites. This process is iterated 
for N times, each time producing a new layer in the bulk 
that has 1/2 D the number of sites as the preceding layer, 
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until only one site is left. The resultant (bulk) tree[35] 
is unitary equivalent to the original (boundary) system, 
and is illustrated in Fig. 1. 



FIG. 1. (Color Online) Illustration of the EHM on 2 3 = 8 
sites. At each iteration, two auxiliary sites (blue) is fed into a 
unitary transform that produces a UV (red) DOF that defines 
a bulk site, and a IR (blue) DOF that becomes the auxiliary 
site for the next iteration. The bulk sites are arranged in a 
tree-like structure (red triangle) with 4 layers, inclusive of the 
last (lowest energy) IR that forms the last ”layer”. 


B. Emergent bulk geometry through boundary 
correlators 

The key motivation behind the EHM approach is to 
uncover the relationship between space-time geometry 
and the quantum entanglement properties of a quantum 
many-body system. The unitary mapping defined by the 
tensor network defines a new direct-product decomposi¬ 
tion of the Hilbert space , and is chosen to make physical 
correlation functions more local in this new basis. To be 
more precise, we assume that the two-point connected 
correlation functions in the bulk always decay exponen¬ 
tially, according to the geodesic distance of certain emer¬ 
gent geometry: 

C f (xi,ti)(x 2 ,t 2 ) — (Oxi( g)o X2 (£2)) — ( 0 Xl (G)) {C X 2 (^2)) 
oc exp ( ^(x 1 ,t 1 ),(x 2 ,t 2 )/ < C) (3) 

This assumption can conversely be used as a definition 
of the distance [30]: 


a — t ^( x i4i)(x 2 ,t 2 ) (A ^ 

d( Xutl ),(x 2 ,t 2 ) - -4 iog- 77 - (4) 

Go 

where Co and £ control the overall offset and scaling re¬ 
spectively. £ can be physically interpreted as the inverse 
mass of the emergent bulk theory, which may depend 
slightly on how we perform the EHM. The logarithmic 
dependence is physically motivated by the observation 
that for a massive system, d( Xljtl ) ? ( X2)t2 ) should recover 
the Euclidean distance in the original system. 

In this work, we shall for simplicity focus on systems 
that are translationally-invariant in space and time, and 


study only correlators with purely space or time intervals, 
i.e. At = £2 — £1 = 0 or Ax = — x\ = 0. 

In the former case with purely spatial interval, all two- 
point connected correlators are bounded above [36] by the 
mutual information 

/ X y = S X + Sy ~ S X y (5) 

where S x and Sxy are the entanglement entropies (EE) of 
a single site and two sites respectively. Roughly speaking, 
the mutual information between two sites measures how 
much the entanglement entropy of two sites will be re¬ 
duced if the correlation between the two sites are known. 
Hence a basis-independent definition of the spatial geom¬ 
etry is given by the mutual information: 

4x = -£ log G (6) 

where Ax = |x — y| and Iq is a reference value for 
the mutual information. One reasonable choice of Jo 
is Jo = 2Klog2, which is the maximal mutual informa¬ 
tion between two sites, each with K internal DOFs (spins, 
bands, etc). This bound is saturated in the (hypotheti¬ 
cal) situation when S xy = 0 but S x = S y = Hlog2, i.e. 
when the DOFs of the two sites are maximally entangled 
with each other but not with those of the other sites. A 
more detailed explanation for the mutual information is 
given in Appendix A 2. 

In the latter case with purely temporal interval, we 
specialize Eq. 4 to 

d ’ = - iAog W) 

where C(r) represents chosen component/s of 
(O x (0)O x (r))*mZfc, r being the imaginary (Wick- 
rotated) time interval r = —iAt. The imaginary 
time direction is preferred over real time as the latter 
typically exhibits oscillatory behavior that makes an 
asymptotic comparison difficult. Further discussion on 
the relationship between the real and imaginary time 
correlators will be deferred to Appendix H. Note that 
unlike the case with spatial intervals, there is no known 
operator that yields the upper bound of correlators 
across temporal intervals. 

Henceforth, we shall use Eq. 6 and 7 as the expressions 
for the distance between two points in the bulk, and com¬ 
pare them with the geodesics of classical geometries. 


III. EXACT HOLOGRAPHIC MAPPING FOR 
FREE LATTICE FERMIONS 

We now specialize the above developments to free lat¬ 
tice fermions, for which the correlators and mutual in¬ 
formation possess nice analytic behavior, at least asymp¬ 
totically. First, we recall the following well-known result 
for the entanglement entropy of free fermions [37-39]: 

Sx = -Tr (Cx log Cx + (I — Cx) log(I - Cx)) (8) 
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where Sx = — Tr (pxlogpx) is the entanglement en¬ 
tropy for the region V, and Cx is the projector (cor¬ 
relator) onto region X. With the help of Eq. 8, it is 
shown in Appendix A 1 that the Mutual Information is 
approximately 


/xy = Sx + S : 


y ^xy 




1 


Cx-y n n n x Cy-x + (x <+ y) 


C y (I-C y ) 


' [^y-x^y-x] 


(9) 


where C x , C y are the single-particle onsite correlators, 
and C x - y is the single-particle propagator between the 
two different sites x and y. This result is completely 
general, and implies that 


4x = (ax log ~ 2 £ax log Tr C y _ x (10) 

in the limit of large spatial separation |x — y|. That 
C x ,C y drops out is hardly surprising, as they each de¬ 
pend only on one site, and have no knowledge about their 
separation. Indeed, most of the information transfer in 
the asymptotic limit is dominated by the single-particle 
propagator. 

We also define the temporal distance via 

d T = ~^ T I 0 ® Tr C(0) (U) 


where a trace of the fermion states have been taken. This 
is the simplest possible basis-independent combination of 
the components of C(r). 

In the next two subsections, we shall introduce pro¬ 
totypical fermionic models as the boundary systems in 
1 + 1 and higher dimensions, and show how their cor¬ 
responding bulk distances and hence geometries can be 
computed via suitable holographic unitary mappings. 


A. EHM for (l+l)-dimensional lattice Dirac 
fermions 

The (l+l)-dimensional lattice Dirac model is among 
the simplest models with a single critical point. In this 
subsection, we will summarize the explicit construction of 
the EHM for this system. Its simplicity allows us to study 
its multitude of entanglement and geometric properties 
analytically with minimal complication. 

The (1 + l)-dim Dirac hamiltonian is a 2-band hamil- 
tonian given by 

Hoiracik ) = Vp[smkcTi + M(m + 1 — COS&)(J 2 ] (12) 

where < 71,02 are the Pauli matrices and vp, the Fermi 
velocity, controls the overall scale of the dispersion. M 
is controls the relative weight between the sin k and 
m + 1 — cos k terms, and will be set to unity here. A 
discussion for generic M will be given in Appendix F. 


When m = 0 or ±2, its gap closes at k = 0 and it be¬ 
comes critical with two crossing bands with linear dis¬ 
persion. To explore or ’’zoom into” the low energy (IR) 
degrees of freedom (DOFs), we utilize a unitary trans¬ 
form that maps states l^ 2 ) on neighboring sites 

into symmetric (low energy) and antisymmetric (high en¬ 
ergy) linear combinations (IV’i 1 ) =b IV^ 2 )- Note that 
the unitary transform does not rotate the spin labels 
si, S 2 , which we shall suppress in the following. In matrix 
form, the unitary transform is written as 

u -+(i - 1 ) (13 > 

The symmetric combination has a Fourier peak at k = 0, 
which is exactly the gapless point of the critical (m = 0) 
Dirac model. The discerning reader will notice that U 12 
is nothing other than the defining expression for the Haar 
transform. Indeed, the construction of the EHM ba¬ 
sis is mathematically identical to performing a wavelet 
decomposition[40-42]. A systematic study of all possi¬ 
ble wavelet descriptions of the EHM will be deferred to 
future work, since for this work we will be primarily con¬ 
cerned about the behavior of the bulk geometries due to 
qualitatively different boundary systems, not the details 
of the wavelet mapping. The transform given by Eq. 13 
possess the virtue of simplicity and, most importantly, 
fixes the archetypal Dirac Hamiltonian, a property we 
shall prove in the next subsection. 

More insight into the EHM can be gleaned in momen¬ 
tum space, where one can directly see how the Hilbert 
space is decomposed into layers with different momen¬ 
tum spectral distributions. Fourier transforming the ac¬ 
tion of Eq. 13 on the single particle states, we obtain 

K) = T, 2 k c ( eik )\^ 2 k) and 1 0k) = T, 2 k D ( eik M 2 k) for 
the auxiliary and bulk states respectively, where |V+) is 
the periodic part of the Bloch state and 

C(e ifc ) = + ( l+e+ (14) 

D(e ik ) = + (1 - e ik ) (15) 

We shall call C, D the IR and UV (low energy and high 
energy) projectors. Physically, they represent the spec¬ 
tral weight projected to the auxiliary and bulk DOFs 
at each iteration. Through these iterations, we ob¬ 
tain successive basis projectors for each bulk layer that 
are increasingly sharply peaked in the IR. To under¬ 
stand this, note that the basis projector of the n th layer 
W n (z) = W n (e lk ) is obtained from n — 1 consecutive IR 
outputs \a) and one final UV output \/3). Hence the first 
bulk layer should contain the DOFs projected from the 
UV projector D(e lk ), while the second layer should con¬ 
tain an IR projector C(e lk ) followed by an UV projector 
D(e 2lk ) that peaks at half the momentum. This rea¬ 
soning generalizes to all the N layers, so the normalized 
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projector for the n th layer is given in momentum space 
by (writing z = e lk ) 


n—1 


i=i 

2 " . 

n 


V2tt 

1 1 — z 


i +Z*- 1 ) 


i i 


(i -z^y 


VzttVz™ 


l-z 


(16) 


W n (e lk ) contains a series of peaks interspersed by valleys 
at e l2Tl k m 1. The dominant peaks occur at k = ±&o ~ 
where the denominator is most singular, as shown in 

Fig. 2, and has magnitude \W n (e lk °)\ = This 

means that as n increases, the spectral weight of the n th 
bulk layer exponentially approach the IR point at k = 0. 
One can further show that the W n ’s form a complete 
an orthonormal basis, i.e. W^ l (e ik )W n (e 2 ' k )dk = 

§\z\~i w m( z ~ 1 ) w n(z)^f = S mn , where the conjuga- 
tion symbol in W* denotes that only the coefficients of 
W(z ), not the argument z, are complex conjugated. In¬ 
deed, that the W n : s are orthonormal with peaks ko ~ 2 _n 
is testimony to the fact that the EHM is a unitary map¬ 
ping that separates the momentum (or energy) scale. 

Note that the auxiliary projector, i.e. projection to 
auxiliary sites with the IR (low energy) half of the DOF, 
is just the orthogonal complement of W n (z ) in Eq. 16: It 
is given by ^== nu c ^ 2 ), comprising IR projectors 
C(z ) only. 



1 2 3 4 5 6 k 


FIG. 2. (Color Online) Plot of the normalized spectral weight 
\W n (e^)\ 2 of the bulk basis for n = 1,2, 3,4. We see that 
as n increases, the dominant spectral peak approaches the 
unique IR point at k = 0 (and its periodic image at k = 2i r) 
exponentially viz. ko ~ —»• 0. Also, it becomes narrower 
since its spectral width also goes like ~ 2 -n . W\ is peaked 
at the highest momentum k = tv, attesting to the fact that it 
contains all the UV DOFs. 

From Fig. 2, we see that the bulk basis from Eq. 16 


become more effective in separating different momentum 
(and hence energy) scales as n increases, since they be¬ 
come more sharply peaked. This property is due to the 
recursive nature of the definition Eq. 16, which also im¬ 
plies that physical quantities, i.e. correlators must scale 
universally with n in the IR limit. However, departures 
from universal scaling may occur at small n (UV regime) 
due to the non-universal high energy characteristics of 
the boundary system. 

Note that the abovementioned unitary EHM transform 
successively ’’zooms into” the low energy DOFs of any 
model with a critical point at k = 0, and not just that of 
the Dirac model. 

Having defined the mapping explicitly, we now write 
down the explicit expression of the bulk correlators. Pre¬ 
viously, we have seen how the bulk distance can be ex¬ 
pressed in terms of the bulk correlators C x ,C x - y and 
C(t), which are the onsite, spatial and temporal propa¬ 
gators respectively. In the free fermion system discussed 
here, the bulk correlators are determined by the bound¬ 
ary two-point correlators 

G q (r) = + e W(ci)-^-i ( 17 ) 

where H( q) is the boundary single particle Hamilto¬ 
nian matrix, (i the chemical potential and T = 1//3 the 
temperature. Explicit expressions for G q (r) as well as 
their resultant mutual information / xy are derived in Ap¬ 
pendix B For cases with particle-hole symmetry, which 
includes the Dirac model. In the zero temperature limit, 
G q (0) reduces to a projector onto the occupied bands be¬ 
low the chemical potential. In the extremely high tem¬ 
perature limit, it becomes nearly the identity operator, 
which just means that almost every state is equally ac¬ 
cessible. In the remainder of the paper, we will always fo¬ 
cus on imaginary time correlator unless otherwise stated, 
since it is still unclear how to define time-drection dis¬ 
tance from the rapidly oscillating real time correlators 
(as discussed further in Appendix H). 

The bulk correlators are most easily expressed as a 
sum in momentum space, since we have already found 
projections to the various bulk layers in terms of the 
spectral weight. Taking the thermodynamic limit where 
the boundary system is infinitely large, i.e. N oo, 
the sum over momenta can be replaced with an integral. 
The bulk correlator between two bulk points (#i,ni,0) 
and (x 2 ,ti 2 ,t) is given by 

C(ni,n 2 , A(2 n x),r) 

*= Y, W nM q ) W n^ iq Y qA{ - 2nx) G q {r) (18) 
q 

where A(2 n x) = 2 n2 x 2 — 2 ni xi is the bulk angular inter¬ 
val and G q {r ) is the (matrix-valued) (1 + l)-dim bound¬ 
ary correlator. The indices ni, n 2 specify the coordinates 
(’layers’) in the emergent ’radial’ momentum-scale direc¬ 
tion in the bulk. There is no translational symmetry 
in this momentum(or energy)-scale direction, unlike the 
original spatial and temporal directions. 
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We see that C(rii, 712 , Ax, r) is just a Fourier transform 
of the boundary correlator G q weighted by the spectral 
contributions W* W n2 of the respective bulk layers. One 
important observation is that the spatial coordinate in 
the exponential factor is A(2 n x), not Ax. This is be¬ 
cause the correlation comes from the holographic projec¬ 
tion of the bulk sites onto the boundary, which depends 
on the angle subtended by the bulk displacement: In the 
simplest case of a circular boundary, the angle subtended 
by Ax at the n th layer is AO = = qq(2 n Ax). The 

bulk correlator still possesses translation invariance, but 
of A(2 n x), the angular interval projected on to the bound¬ 
ary , not of the bulk interval Ax itself. Mathematically, 
we find that the 2 n rescaling of angular distance is also re¬ 
quired for the orthogonality of the basis {W n (q)e l2 q ' x }, 
x = 1 , 2 ,..., 


B. Transformation of the Hamiltonian under the 
EHM and its fixed points 


The EHM is an exact version of renormalization group 
(RG) transformation. In each step, the DOFs of the sys¬ 
tem are split into the high energy (UV) and the low en¬ 
ergy (IR) parts, with the procedure iterated on the low 
energy part. If we write the Hamiltonian in the new ba¬ 
sis after an EHM step and ignore the coupling between 
the IR and UV degrees of freedom, we can write down 
a ’’low energy effective Hamiltonian” of the IR states. 
This resembles the renormalization group flow of the ef¬ 
fective Hamiltonian in ordinary RG. For a given choice 
of the EHM tensor network, there are certain bound¬ 
ary systems for which the IR Hamiltonian is at an RG 
fixed point. In the following, we will show that the mass¬ 
less Dirac Hamiltonians are RG fixed points of the EHM 
transformation we defined earlier. 

We start by writing down the effective IR Hamiltonian 
in the momentum basis. In D = 1 spatial dimensions, the 
EHM for each iteration is given by the change of basis in 
Eqs. 14 and 15, so the single particle Hamiltonian matrix 
h n of the n th layer is transformed according to 


h n+ l( e i*/2) = 


yt f h n (e ik/2 ) 


11 


(19) 

with the two components of the matrix represent¬ 
ing the IR and UV DOFs. Here V(w) = 
i ( C(w) D(w) \ _ i f 1 + w 1 — w 
U2 y C(—w) D(—w) J 2^ 1 — w 1 + w 

w = e lk ! 2 . To obtain the Hamiltonian at the (n + l) th 
layer, one projects onto the upper left or IR component 
of the RHS. With D > 1 spatial dimensions, V will be 
given by V(w\) 0 ... 0 V(wd). If the Hamiltonian is to 
remain invariant under the RG, the 11 (IR) component 


^, where 


in Eq. 19 gives 2 h n+1 (w) = 2A h n (w 2 ) or, in detail: 

2A h n {w 2 ) 

= + h n (-w)C(-w)C*(-w^) 

= (h n (w) + h n (-w )) + W + W 1 (h n (w) - h n (-w)) 

( 20 ) 

where A is a constant scale factor for each EHM step. 
Upon setting w = 1, we obtain 

Aft(l) = h( 1) (21) 

which implies that A = 1 unless h( 1) = 0, i.e. that the 
rescaling A for each step can be nontrivial (A / 1) only 
if the Hamiltonian is gapless [43] at the IR point k = 0 
or w = 1. In other words, only gapless Hamiltonians can 
have nontrivial scale invariance, as is expected. This has 
very important implications in spatial dimensions D > 1, 
since it implies that if the Hamiltonian contribution h 
does not depend explicitly on kj , the Hamiltonian will 
not be scale invariant by any EHM transformation in the 
j th dimension. 

One can derive solutions of Eq. 20 by comparing terms 
power by power. For instance, the second-order terms 
in h(w 2 ) on the LHS force h(w) to be a linear function 
of w that is either symmetric or antisymmetric under 
w «+ re -1 , i.e. a function of w-\-w _1 or w—w -1 . Hence we 
find the two linearly independent solutions to Eq. 20 to 
be h(z) = = sin k and h(z) = 2 ~^~ 1 = 

both with the EHM rescaling A = ^. They are both 
gapless at k = 0, as they should be, and can be combined 
to form the massless Dirac Hamiltonian Eq. 12 in 1 + 1 
dimensions: 

Hoiradk ) = up [sin k(Ji + M( 1 — cos k)<j 2 \ (22) 

where cq are the Pauli matrices and M controls the rel¬ 
ative weight of the two terms. 

Note that the two terms sin k and 1 — cos k do not 
have the same scaling dimension if one takes the contin¬ 
uum limit sin k ~ k, 1 — cos k ~ k 2 /2 in ordinary RG. 
This illustrates the distinction of real space EHM trans¬ 
formation from simple momentum rescaling, due to the 
nontrivial influence of lattice regularization that replaces 
functions in k -space with periodic trigonometrical func¬ 
tions. 


IV. ANALYTIC RESULTS FOR 
(1 + l)-DIMENSIONAL BOUNDARY SYSTEMS 

In Ref. 30, the behaviors of correlation functions 
and their associated bulk distances were studied numer¬ 
ically. In the following, we will obtain analytic results of 
the asymptotic behavior of bulk correlators for (1 + 1)- 
dim translationally-invariant boundary systems, which in 
turn determine the asymptotic large scale behavior of the 
bulk geometry we define. We will compare them with the 
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geodesic distances between analogous points in candidate 
classical geometries. The higher-dimensional extensions 
of these results will be discussed in the next section. 


A. General setup 

We shall consider four distinct physical scenarios, all at 
zero chemical potential, with the representative Hamilto¬ 
nian for the first three cases taken to be the Dirac Hamil¬ 
tonian given in Eq. 12. The results obtained should also 
be valid for more generic Hamiltonians, since the qual¬ 
itative bulk geometry properties remain robust as long 
as the long distance behavior of correlators remain the 
same. In approximately increasing levels of sophistica¬ 
tion, the four scenarios are: 

1 . Critical boundary Dirac model at T = 0, corre¬ 
sponding to a bulk AdS (Anti-de-Sitter) geometry. 

2. Massive boundary Dirac model at T = 0, corre¬ 
sponding to a “confined geometry” with an IR ter¬ 
mination surface. 

3. Critical boundary Dirac model at T ^ 0, corre¬ 
sponding to a bulk BTZ (Banados, Teitelboim, and 
Zanelli) black hole geometry. 

4. Critical boundary model with nonlinear dispersion 
at T 7 ^ 0, corresponding to a bulk Lifshifz black 
hole geometry. 

The first two cases were already explored numerically 
in Ref. 30, with results in excellent agreement with our 
analytical results below. 

As previously explained, the fundamental quantity to 
be calculated is the bulk correlator C(ni, 712 , A(2 n x), r) 
given in Eq. 18. In the thermodynamic limit L —>• 00 , all 
momentum sums can be replaced by integrals: 

C(n 1 ,n 2 , A(2 n x),r) 

= Y, W nM) W nMe iqA ^ x) G q {r) 
q 

= <f -W* ni {z^)W n2 {z)z^ 2 ^G z {r) 

J\z\ = l z 

(23) 

where, as before, the conjugation symbol * indicates that 
only the coefficients of the polynomial W(z) are complex 
conjugated. It is insightful to analytically continue the 
momentum q into the complex z = e iq plane, where the 
decay properties of the correlators can be directly read 
from the properties of the complex poles and branch cuts. 

For comparison with the geodesic distances, we shall 
specialize to 2-point correlators of the following three di¬ 
rections in the (2 + l)-dim bulk: 

• Equal time, same layer “angular” correlator 
C n {Ax) = C(n, n, Ax, 0) 

= <f -W:{z~ 1 )W n {z)z 2nAx G z { 0) (24) 

J\z 1 = 1 * 


• Equal time, different layer “radial” correlator 

C(n 1 ,n 2 ) = C(n i,n 2 ,0,0) 

= <f ^WZiz-^Wn^GM (25) 
J\z\=l z 

• Same site imaginary-time correlator 

C n {T) = C(n, n, 0, r) = <f -W*(z- % )W n {z)G, z {t) 

J\z\=l z 

(26) 

Here, we have assumed translational invariance in x, 
which is necessary for defining the correlator in terms 
of a Fourier integral in the angular direction. 

Next we specify the boundary Hamiltonian. We shall 
use the Dirac Hamiltonian 

HDirac(k) = VF[sinkai + M(m + 1 — COS&)<72] 

from Eq. 12 for cases (1) to (3). For the sake of concise¬ 
ness in the already sundry results, we shall henceforth 
set vp and M to unity unless otherwise stated, and con¬ 
sider only cases at zero chemical potential fi. Indeed, xp, 
which couples to r under imaginary time evolution e~ Hr 
merely leads to a trivial rescaling r vpr in the results. 
The value of M does not affect the leading asymptotic 
behavior of the correlators in general, and its study is rel¬ 
egated to Appendix F. For case (4), we shall simply base 
our calculations on the non-linear dispersion relation 

E k = k\ 

since we will be primarily interested in the effect of set- 
ting 7 ^ 1 . 

We next elaborate on the complex analytic structure 
of the Dirac Hamiltonian and correlator. The positions 
of the complex singularities play a crucial role in deter¬ 
mining the asymptotic decay behavior of the correlators, 
typically with power-law decay when all singularities he 
on the unit circle and exponential decay otherwise. The 
correlator in the spatial ” angular” direction, in particu¬ 
lar, is a Fourier transform for which there exist results 
that relate the decay of Fourier coefficients with the lo¬ 
cation of singularities. For a meromorphic function f(z), 
the Fourier coefficients fi = §^ =1 ^ff(z)z l decay like 

fi~r (1+B) \z 0 \ l (27) 

for |zo | < 1, / 1, with zo the branch point of f(z) clos¬ 

est to the unit circle and B is its corresponding branching 
number: 

/A \ B 

f(z 0 + A z) ~ f(z 0 ) + (-^-J ( 28 ) 

for z near zq. Note that B cannot be a non-negative 
integer, since otherwise the Riemann surface will not be 
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ramified or even divergent at zq. This result from com¬ 
plex analysis has been heavily used in a variety of physi¬ 
cal problems, from the decay properties of states in con¬ 
densed matter to the bulk properties of statistical net¬ 
works. Proofs, together with its physical applications, 
can be found in Refs. 44-47 and especially 48. 

In our correlators of interest, f(z ) takes the explicit 
forms h z or h z /E z from the expressions to follow. h z is 
the Dirac Hamiltonian with M = l : vp = 1: 


h z 


0 

—i(z — (1 + m)) 


i(i-(l + m+j 


(29) 


with eigenenergies 


E z = 1 + (m + l) 2 - (m + 1) ( z + - 


m —0 


» —1{Z 2 — 2 2 ) 


(30) 


Due to the square root, 2 is an analytic function on 
a 2-sheeted Riemann surface with ramification (branch 
points) at z = 00 , m + 1, ^-j- and 0. When m = 0, the 
two points z = (1 + m) ±:L coincide and annihilate, leaving 
a single branch cut from 0 to 00 . These branch points 
also appear in the flattened Hamiltonian that appears 
in the correlators Eqs. B 6 , B7 and B 8 : 


/ 


hz 

E z 


m+1 


2 —(ra+1) 


m+1 


z— (m+1) 

m+1 


m—>-0 


0 


+ 0 / 


(31) 


When m = 0 at criticality, takes a particularly simple 
form that does not contain any nonzero pole in the unit 
circle. This still holds true for much more generic critical 
systems, being a necessary condition for power-law decay 
as required by Eq. 27. Note that all the nontrivial branch 
points in the integrand of the bulk correlator must come 
from G(z), since the wavelet basis functions W n (z) or 
W n (z~*) are polynomials in z or 

Interestingly, this simple form of |+ admit a real 

space correlator G(Ax) oc f z Ax ^E with a nontrivial 
phase winding. Indeed, a short calculation reveals that 

G(Ax) oc (aI+i/ 2)] ’ exac % agrees with nu¬ 

merical results for different models with a Dirac point [49]. 

Although we have only explicitly studied the complex 
topology of the Dirac model, the important point is that 
more generic models, i.e. multi-band models with arbi¬ 
trary dispersion also have analogous topologies that lead 
to similar asymptotic correlator behavior. This will be 
further elaborated in the last part of Appendix F. 

Table I contains a summary of the asymptotic behav¬ 
ior of the mutual information, correlators, entanglement 
entropy and bulk geometry parameters for the different 
boundary systems. More details about each case are dis¬ 
cussed in the remainder of this section. 


B. Critical boundary Dirac model vs. bulk AdS 
space 


For a start, let us consider the massless (1 + l)-dim 
massless Dirac model as the boundary theory. We show 
that the correlators in all three directions (Eqs. 24 to 
26) all suggest a bulk geometry of a ( 2 +l)-dim Anti de- 
Sitter space (further detailed in Appendix I) given by the 
metric (Eq. 12): 


p 2 j n 2 


L 2 

4t t 2 R 2 


1 4 - El 
1 E R 2 


(32) 


where L is the boundary system size and r is scaled such 
that the metric is 0 ( 2 )-invariant at the critical boundary 
where p = -E r l n the following, we shall only con¬ 
sider asymptotically large angular and temporal intervals 
with A x,r -+ 00 , and/or not necessarily large radial in¬ 
tervals between layers 1 and n. Since we have already 
taken the thermodynamic limit L -+ 00 , Ax < L can 
still be satisfied for arbitrarily large Ax. 


1 . Spatial directions 


As is detailed in Appendix B, the spatial decay of the 
Mutual Information is given by (see Eq. (B4)): 


4(M 2 + H 2 ) 

1-4 A 2 


(33) 


where |^|, \v\ are the off-diagonal (unequal spin) part of 
the propagator C x - y , and A is the off-diagonal part of the 
onsite propagator C x . From Eq. 29, they are explicitly: 


A = -iC n {Ax = 0) = i l —W*(z- l )W n {z)L~z 
+ 1=1 * 


(34) 

which, being onsite, does not depend on the displacement 
between x and y. For the angular direction where x,y 
are on the same layer, 


u,v = -iC n (Ax) = —i d) -—W*(z~ 1 )W n (z)z T ^z 2nAx 

j\z \=1 2 

( 35 ) 

while for the radial direction where n 1 7 ^ 712 but Ax = 0 , 


u,v= -iC{n u n 2 ) =-if —Wl^z- x )W n ^z)z^ 

+|=i 2 


(36) 

In the above equations, the integrands do not contain 
poles. However, the integrals are nonzero due to the 
branch cut from z = 0 to z = 00 from the square root 
factor. They can be evaluated by standard deformations 
of the contour, as demonstrated in detail in Appendix C. 

After some computation, we obtain for the angular di¬ 
rection 


M+1 


1 1 

167T ( Ax ) 3 


(37) 
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m — T — 0 (Case 1) 

m/ 0 ,T = 0 (Case 2 ) 

T ^ 0, m — 0 (Cases 3,4) 

s x 

log 4, Eq. 40 

oc ^ 0, Eq. 59 

log 4 

In (Ax) 

~ Eq. 41 

- e a, 2 . E q- 58 

~ e ~ 2nT2rlAx for 7 = 1 ; and in general 

e -[ 2 (,T) 1 /l ( 2 »Ax)]si„^ ) Eq . 8 g 

7(711,712) 

~ 2 -An , Eq. 42 

unexplored 

unexplored 

Tr C„(t) 

o3n „ , „ 

~ ’ Eq - 49 

- e -mT 

Tr [C n (r) + Cn(/3 — t)]t =0 for 2 n <C 27r/3; 

bounded below by e 2 ( 27 - 1 )™ ? gq S go,89 for 2 n > 2^/3 

R 

1 

3 

n.a. 

1 for 7=1 

R 

1 

n.a. 

1 for 2 n 27r/3; n.a. otherwise 

ix. 

R t 

2 

3 

n.a. 

2 

3 

b 

n.a. 

n.a. 

L^qT for 7 = 1 ; oc T 1 ^ 1 for 7 > 1 

Re 

1 

1 / i /TT 3 

2 V ^ V i° ) 

n.a. 

not uniquely determined 

Rt 

( 2 i> F 7 T 1/3 ) -1 , Eq. 51 

n.a. 

\f\ for 7 = 1 

y 7T 1 

Pn 

L 

2 n + 1 7T 

n.a. 

5 ^ for 2 " < 

b ( x + sf7) for 2 " > 


TABLE I. Summary of results for various physical quantities in all of the D — 1 scenarios considered, where L is the size of 
the boundary system. Case (1) ( m = T — 0) from Sect. IV A is fitted onto the AdS metric with radius being either R or R r 
(which are numerically close) depending on whether the fit is done for the spatial direction or imaginary time direction. The 
massive (m ^ 0) case (2) is not fitted with any classical geometry. For cases (3) and (4) with temperature 0, 7 controls the 
dispersion via E q = q 1 , and results are given for the BTZ case (7 = 1) and the general Lifshitz case (7 > 1), when applicable. 
The energy scale of T divides the bulk into two regions demarcated by 2 n = 2tt/3, each with different qualitative properties. 


As elaborated in Appendix C 1, such a power-law decay 
of the single-particle propagator is generally expected in 
the presence of a branch cut. Physically, it is a signature 
of criticality, with a power of 3 instead of 1 due to the ad¬ 
ditional ’destructive interference’ from the antisymmetric 
combinations of adjacent sites in the Haar wavelet basis. 
There is a striking absence of the layer index n in Eq. 
37, which reflects the scale-invariance of the boundary 
theory. Eq. 37 also holds for general values of M in the 
Dirac Model Eq. 12, where M controls the ratio between 
the quadratically dispersive 1 — cos k and linearly disper¬ 
sive sin k terms near the IR point. While M can affect 
the details of the branch cut, it cannot change the decay 
exponent, as is shown in Appendix F. 

For the radial direction with n\ — 1 and 712 = n, Ap¬ 
pendix C 2 also tells us that 


71 — i 

(38) 

Hence we have exponential decay of the single-particle 
propagator in the radial direction, which is consistent 



with scale invariance [50]. 


Strictly speaking, the mutual information across the 
radial direction involves both A\ and A n , the unequal 
spins onsite propagators A at layers 1 and n, and a more 
general (and complicated) version of Eq. (B4) should be 
used. However, A n 0 rapidly as n increases, effec¬ 
tively leading to no asymptotic correction. Mathemati¬ 
cally, this is because as n increases, the peaks of W n (e iq ) 
approaches a delta function at qo = > 0 which gets 

exponentially closer to the IR point, where contributions 
to the integral are penalized by the momentum correlator 
h q /E q . Explicitly, 
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Kl = 


in the large n limit. Physically, this means that unequal 
spins become totally decoupled in the IR regime. Ac¬ 
cording to Eq. B5, the sites in the IR layers are hence 
maximally entangled with the rest of the bulk: 

Sx -» -4^ log i = log 4 (40) 

This maximal entanglement in the IR also exists in 
generic critical systems, since the IR DOFs become 
harder and harder to isolate unless an energy scale (i.e. 
mass) exists. 

Putting it all together, the mutual information behaves 
like 


~lf \W n (e iq )\ 2 sgn(q)e iq / 2 dq 
\ f° q \W n (e iq )\ 2 ^dq 

i r qo 

- A \W n (e iq °)\ 2 / \q\dq 

4 J — cm 


on —1 
z 2 
—3-00 
7 T 6 

2 1 

— 7 \— —^ 0 
7T 2 n 


(39) 


— log/ n (Ax) ~ 61og Ax + log(327r 2 ) (41) 

for angular intervals and 

— log 7(1, n) ~ (n — 1) log 2 + small const. ~ An log 2 

(42) 

for radial intervals. These asymptotic behaviors are in 
excellent agreement with those of the geodesic distances 
on AdS space, if one uses the proposed correspondence 

dmin j 

given by = — log in Eq. 6 . The parameters £ 
and Jo respectively set the scales of bulk distance and 
mutual information, and are related in a precise way dis¬ 
cussed later. In principle, £ can be different in different 
independent directions. 

We first study the correspondence in the angular di¬ 
rection. The geodesic distance between two equal-time 
points (p, 61 ) and (p, 62 ) with angular interval A 6 = 
\02 — 0i\ = is given by Eq. 14: 

dt: = # cosh - 1 (l + ^ sin 2 ~ 2R log ^ (43) 

where R is the AdS radius that determines the length 
scale below which we expect significant deviations from 
logarithmic behavior. Comparing Eqs. 41 and 43, we see 
that 


R 
f e 


= 3 


(44) 


and 




(45) 


We see that R decreases weakly with increasing Iq. For 
the Dirac model we take Iq = 2 log 4, the theoretical max¬ 
imal mutual information mentioned below Eq. 6 . With 
it, we obtain R = 0.3233 and £<9 = 0.1078, which is in 
excellent agreement with the numerical values obtained 
in Ref. 30. 

For the radial direction, we compare Eq. 42 with the 
AdS geodesic distance (Eq. 15): 

d™ n =R\A(logp)\ 

= R(log 2 )|An| (46) 

Here we have taken p = p n = 2 ^^’ as squired by the 
scale-invariance of radius p and the circumference ^ at 
any layer n. We easily obtain 


i P = R (47) 

which means that the AdS radius is nothing but the ra¬ 
dial length scale £ p of radial geodesics. This is also in 
agreement with the numerical results in Ref. 30. 

It should be noted that the ratio i?/£ is different for the 
angular and radial directions, which is a manifestation of 
the fact that the mapping does not preserve the entire 
conformal symmetry of AdS space (since conformal sym¬ 
metry is only emergent in long wavelength limit and does 
not exist rigorously in a lattice model). The bulk theory 
has scale invariance and a reduced translation symmetry 
with unit cell size varying with the layer index n, but 
the correlations are anisotropic between radial and angle 
directions. 


2. Imaginary time direction 


As previously postulated by Eq. 7, we can deduce a 
classical bulk geometry from the decay properties of the 
imaginary time correlator C n {r) given by 

Cn(r) = i dq\W n (e iq )\ 2 e~ TE ' (l - (48) 


When r is large, the value of C n (r) arises from com¬ 
peting contributions from the IR regime and the mo¬ 
mentum scale set by the layer index n. While most of 
the spectral weight of W n {e iq ) is concentrated around 
the momentum go = , the exponential factor exponen¬ 

tially suppresses contributions above the IR regime set by 
E q = vpq ~ t -1 . Hence we expect C n {r) to increase as n 
goes deeper into the IR. Indeed, this is exactly contained 
in the analytical result Eq. E 6 derived in Appendix E: 


Tr C n (r) - 


1 

87 r 



( 87 r 2 ) 2 \vppr 


3 


(49) 
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where p = 7 ^y is the bulk radial AdS coordinate of layer 
n. Eq. 49 also holds for generic M in the Dirac model 
(Eq. 12), except for the degenerate case where M = oo 
and there is no linearly dispersive term sin ka i. 

Comparing the logarithm of Eq. 49 with the AdS time¬ 
like geodesic d T ~ 2Rlog from Eq. 16, we obtain 


R _ 3 

and 

1 ~ 0.3414 

T 2 U^TT 1 / 3 Vp 


(50) 


(51) 


These results are valid in the regime where p y> R r and 
R r « r <C T, i.e. when the geodesics do not come close 
to circumnavigating the AdS space. The AdS radius 
R r obtained here is slightly different from R obtained 
through spatial geodesics in Eq. 45, with the latter con¬ 
taining a weak dependence on the reference mutual in¬ 
formation I~ 1/6 . 

We note that ^ = | obtained from the imaginary 
time correlator is exactly half of that of R = 3 obtained 
from the (spatial) mutual information. This is not an 
inconsistency, but a manifestation of the anisotropy in 
the definitions: the mutual information is asymptotically 
quadratic, not linear, in the two-point correlator, and its 
associated geodesic distance must be doubled. 


C. Boundary Dirac model at nonzero mass m, T = 0 

This is our first and simplest example of a non-critical 
system. A nonzero mass scale m introduces a pole in the 
correlator, which leads to the exponential decay of cor¬ 
relation functions. As we will discuss in this subsection, 
the dual geometry has a spatial ’’termination surface” 
which makes the spatial geometry topologically different 
from that of critical fermions. The temporal geodesics, 
however, exhibit no unusual behavior at the termination 
surface, so the surface is not a black hole horizon but a 
purely spatial cutoff. 


1. Angular direction 

As previously discussed, the geodesics in the spa¬ 
tial angular direction depend on the mutual information 
I xy = w here u and v are the unequal-spin prop¬ 

agators between sites separated by an angular distance 
of Ax within layer n , and A is the unequal-spin onsite 
propagator. 

We first look at how u and v differ from those of the 
critical case. They are given by the off-diagonal compo¬ 
nents of 


—i 


dz 

z 


W*(z~ 1 )W n (z)z 2nAx 


h z 

E~z 


where ^ is the 2x2 matrix given by Eq. 29. For the 

Dirac model with m/0, ^ and hence the integrand of 
Eq. 52 now has square root branch points away from the 
origin, namely at zq = 1 Am and zq = pUy. Unlike in the 
critical case, the singularity zq = pUy now satisfies \zq\ < 
1 and by Eq. 27 determines the asymptotic exponential 
correlator decay with f(z) = W*(z~ 1 )W n (z) jg-. The 
branching numbers B (Eq. 28) at zo are given by B = =b^ 
for u and u, which shall thus decay (for 2 n Ax y> 1 and 
small m) like 


(2 n Ax)- 1=Fl /2 ^ e -m2 n Ax 

U,V ~ (1 + m) 2nA * ~ (2”Aa;) 1±1 / 2 (53) 

We see that v dominates with a factor of 2 n Ax, and is 
itself exponentially decaying with a subleading power-law 
decay. Hence — log I xy will define a rescaled Euclidean 
distance with length scale given by ^. For a Dirac model 
with a generic value of M in Eq. 12, zo = yyyp should 
be replaced by the pole 




M( 1 + m) — >/l + 2 m M 2 + m 2 M 2 

Mi 1 


(54) 


whose corresponding effective mass is given by \zq\~ 1 — 1 . 

As contrasted with the critical case, the mutual infor¬ 
mation in the massive case also depends nontrivially on 
A , the unequal-spin onsite propagator, as we go deep into 
the IR regime below the energy scale of layer n. Recall 
that at layer n, A is given by 


• r ir roffdiag 

A n = Y 


-dq 


(55) 


For large n, the spectral weight is mainly concentrated 
around qo = , where 

fooffdiag sin q dz i(m + 1 — cos q) 

Eq Eq 

i(m + %r) 
ym 2 + q 2 

* ■ 0 + irf 1 " 3 ) < 56> 

where the sing term had been dropped because it is odd. 
Hence Eq. 55 can be rewritten as 

= \ /* 0 + d " 

~ 5 - 4^1 w ^£y 

1 2 n ql 

2 8m 2 7r 3 3 

_ 1 11 

“ 2 


ld=i 


(52) 


3 m 2 4 n 


(57) 
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This leads to an overall enhancement of x _l A2 ~ 4 n in 
the mutual information, which thus behaves like 

— log/ n (Ax) ~ 2 n+1 raAx + log Ax — n log 2 (58) 


Its leading contribution is proportional to 2 n Ax, which 
leads to a key difference of its dual geometry from that 
of the critical fermion. Since the corresponding geodesic 
distance d/\ x = %AxIn(Ax) is linear in Ax, i.e. like an Eu¬ 
clidean distance, we can obtain the circumference a n of 
layer n simply by taking Ax = 2 Ar_n , the number of sites 
in the whole layer. More rigorously, the circumference a n 
should be measured by first taking equally spaced points 
with 2 n sites between them, and then taking the N —» oo 
limit before taking the N' —» oo limit. The circumference 
will be given by the product of the geodesic distance be¬ 
tween neighboring points 2 m^Ax • 2 n+Ar and the number 
of points 2 N ~ n ~ N \ i.e. a n = 2 m^Ax * % N which is a finite 
portion of the boundary circumference 2 N . 

The fact that a n is finite in the large n (infrared or IR) 
limit tells us that when n oo, we are not approaching 
the center of a hyperbolic disk, but rather approaching 
a surface with finite area which acts as a “termination 
surface” of the space. This is illustrated in Fig. 3. The 
IR surface shrinks with decreasing mass m. By compar¬ 
ison, in the dual geometry corresponding to the critical 
fermion, the distance between points Ax sites apart does 
not depend on the layer index, so that the circumference 
OL n oc 2 ~ n decays exponentially in the IR limit. 

It is important to note that this surface is not a black- 
hole horizon, which distinguishes it from the case of 
nonzero temperature state we will discuss in subsequent 
sections. One evidence for this conclusion is that each 
site carries a vanishing entanglement entropy S x in the 
IR limit n oo, since 


3 m 2 4 n 
log 16 n 
3 m 2 4 n 


1 1 

3 m 2 4 n 


log 


1 


3 m 2 4 n 



which is obtained by substituting Eq. 57 into Eq. B5. 
Since the entropy decays exponentially at large n, the IR 
sites actually for direct product states unentangled with 
one another. This is consistent with the physical picture 
that mass is renormalized to exponentially larger values 
in the IR limit (with respect to the kinetic energy scale), 
which forces the IR ground state to be simply a direct 
product of the single-site ground states with a large mass 
term. 


2. Imaginary time direction 

The correlator in the imaginary time direction 
can be obtained pretty straightforwardly. Substitut¬ 
ing the expression for the energy dispersion E q = 

yj sin 2 q + (1 + m — cos q) 2 « m + (^pr) q 2 into Eq. 48, 


we obtain 


Tr C n (j ) = e~ mT f” dq\W n {e ic >)\ 2 e- T2 A 

J —TV 


oc e 


(60) 


which is an exponentially decaying term multiplied by a 
nonuniversal Gaussian Integral. Hence — logTr C n (r) ~ 
mr defines an Euclidean bulk geodesic in the imaginary 
time direction. Together with the previous results on 
spatial geometry, we see that although there is a spatial 
termination surface in the IR limit, the temporal direc¬ 
tion still extends as usual. This is consistent with our 
earlier statement that the IR termination surface is not 
a black hole horizon, since the time direction is not in¬ 
finitely redshifted. Topologically, the space-time we ob¬ 
tain for the massive Dirac model is IZ x A4, where 7Z is 
the line in time direction, and M is a spatial annulus. 


Critical system 


With mass/temp scale 




Geodesic 


Blue - Approx. AdS 
Red - Near Horizon 


FIG. 3. (Color Online) The spatial bulk geometries of crit¬ 
ical (Left) and non-critical (Right) systems and examples of 
their geodesics. In the AdS bulk (Left) corresponding to a 
critical boundary system, geodesics extend inwards towards 
the center, where there is far less ’space’. With an energy 
scale introduced by nonzero mass or temperature, the bulk 
geodesics develops a spatial ’termination’ of space (Right). 
Geodesics wrap around a ’prohibited’ region of radius oc m 
or T, acquiring an Euclidean character is that regime. As 
explained further in the main text, only the case at nonzero 
temperature corresponds to a true horizon. 


D. Critical linear boundary Dirac model at 
nonzero temperature T vs. bulk BTZ black hole 


Now we study the effect of nonzero temperature in the 
massless Dirac model, the bulk geometry dual to which 
will be interpreted as a BTZ (Bahados, Teitelboim and 
Zanelli) black hole [51]. The BTZ black hole is a black 
hole solution for (2 + l)-dim gravity with a negative cos¬ 
mological constant, whose metric is 


ds 2 


V(p) 

V(L/( 2 n)) 


dr 2 + 


dp 2 

W) 


+ p 2 d0 2 


(61) 
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2_i 2 

where V(p) = p R 2 is the Lapse function, b is the 
horizon radius and R is an overall length scale. We 

have rescaled r by a factor of [V (^)] so that 
ds 2 —> dr 2 + p 2 d 0 2 possesses 0(2) symmmetry at the 
boundary where 2irp = L. The rescaling makes r the 
same imaginary time variable as that in the boundary 
system (the canonical conjugate coordinate to the bound¬ 
ary Hamiltonian). 

Most notably, we find that the layers deep in the IR 
accumulate at the bulk horizon radius p —b from above, 
where b is proportional to the temperature T = yd -1 . 
This result is rigorously accurate in the m ~ 0(1) range 
of the Dirac model given by Eq. 12, where the critical 
dispersion is essentially linear. 


1. Angular direction 


Like a nonzero mass, a nonzero temperature also intro¬ 
duces an energy scale into the system. This energy scale 
is manifested as an imaginary gap — log\zo\ >0 where zo 
is the singularity of the momentum-space correlator zq 
closest to the unit circle. In this case, the singularities 
originate from the tanh term [52] given in Eq. B8. 

Above the energy scale of T, all correlators and hence 
the mutual information do not feel the effect of ther¬ 
mal excitations, and define an approximate AdS geome¬ 
try just as in the T = 0 critical case. 

As one goes below the energy scale of T, Eq. 27 
states that the unequal-spin propagators u , v decay expo¬ 
nentially as \z $\ 2 Ax , with a subleading power-law term 
(2 n Ax)~^ B \ B and \zq\ can be found as follows. 

Expanding cosh , the denominator of the singular 
term, about zq, we find that 


n , P(z-z 0 ) . (3E z dE z 

cosh ——- = 0 H---smh — -—— oc z — zq (62) 

Zi Z Z CLZ 


so the branching number B = —1. This value of 
branching number holds universally for generic systems 
at nonzero temperature [53], since Eq. 62 does not de¬ 
pend on the form of E z . Therefore, there is rigorously no 
subleading power-law term in the decay of correlators in 
the nonzero temperature case. 

We now proceed to find \zo\. tanh is singular when 
cosh = 0, i.e. iir = f3E q = —ifd(z x ! 2 — z~ x ! 2 ) (see 

Eq. 30). This is satisfied by z + i — 2 = ^ 27r (^ +1 ) ^ ? 
where / G Z. Hence the poles occur at 

±_1 , 7T 2 (2l + l) 2 ^ tt( 1 + 2l)y/4p + (2 1 + l) 2 7T 2 

z ° - + 2/3 2 2/3 2 

(63) 

in the limit of small T = ^. The pole with Z = 0 is closest 
to the unit circle, and we thus have the asymptotic decay 

u ~ v ~ \zq( 1 = 0)\ 2nAx (1 — 7 T T) 2nAx e - ?rT2r ' A;E 

(64) 


so that 


o p -nT2 n+1 Ax 

~ 1-4A2 (65) 

Since A n —>> | after the first few n, as explained in the 
subsection on the zero temperature critical case (1), the 

mutual information behaves like (recalling that 2 n Ax = 
LAO \ 

2tt ) 

- log I n (Ax) - 2'kT2 u A x = LTA6 (66) 


with no logarithmic subleading term. This asymptotic 
form for I n (Ax) is, to leading order, the same as that 
of Eq. (58) for the massive zero temperature case, if we 
replace 2m by 2 i tT. Following along the same lines as 
the previous section, we conclude that the circumference 
of each circle at layer n is asymptotically a n 2 ttT • 2 N . 
Thus there is a termination surface with this circumfer¬ 
ence (1-dim area) in the IR (low energy) limit. However, 
as we will verify by the single-site entropy and imaginary 
time direction distance later in this section, this surface 
is not just a termination surface of space, but a black 
hole horizon. Before discussing that, we shall first make 
a detailed comparison of the angular direction distance 
defined by Eq. 66 with the angular geodesic distance 
given by a BTZ black hole metric. 

We can obtain a precise relationship between the tem¬ 
perature T and the black hole radius b. Since Eq. 66 
holds in the IR regime, we compare it with the BTZ 
geodesic distance[54] d™™ = 2Rsinh -1 [£ sinh ^-\ from 
Eq. 23 of Ref. 30 in the near horizon limit 0 < p — b <^b. 


i r / A \ ^ 2iZ _i 

— log I n (Ax) = — = -—smh 

€o & 


which implies that 


p . , bAO 

- smh —— 
b 2 R 


Ell 0g (P p bAe/2R') 


2 R , p 2 n+1 bnAx 

= — log £ + ——— 

ie b £L 

^ 27r&(2 n Ax) 

* eT 


(67) 


b = Li s T 


( 68 ) 


Here £# is a yet-undertermined length scale. Our bulk 
geometry has agreed remarkably well with that of an ac¬ 
tual BTZ black hole, with the LHS and RHS of Eq. 67 
agreeing on not just the leading term linear in Ax, but 
also the vanishing logarithmic subleading term. Further 
discussions of the near-horizon geometry can be found in 
Appendix 12. 

£<9 can be determined by compairing Eq. (68) with the 
relation between T and b in classical gravity. The require¬ 
ment that the geometry is smooth in imaginary time at 
p = b, i.e. without a conical singularity, requires[55, 56] 


T = 



L 


b 2 i tR 


-l 


b 

~RL 


(69) 
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More details about this formula are given in Appendix 
12. Comparing Eqs. (68) and (69) we obtain 

R = Ze (70) 

so that £<9 and R are actually one and the same length 
scale parameter. 


As suggested in Fig. 4, the first term gives a curve with 
geodesic distance qualitatively similar to that outside a 
BTZ black hole where p^> b: 


d T = R cosh 1 



7TT 

T 


(75) 


which is derived in the Appendix of Ref. 30. 


2. Imaginary time direction 

The nonzero temperature features prominently in the 
imaginary time correlator because a finite /3 corresponds 
to a finite periodicity in imaginary time. We shall show 
that the correlator Tr C n {r) defines a bulk geometry that 
is qualitatively similar to that of BTZ black hole, and 
deduce the effective radius p = p n of layer n by looking 
at the maximal value of — logTr C n (r) achieved at half¬ 
period r = |. 

The quantity to be computed is 

Tr C„(r) = i J* dqlWnie^Tr G q {r) (71) 

where, from Eq. 17, 



FIG. 4. —logTr C n {j) according to Eq. 74. Plotted are the 
n = 1, 2 curves for T = 0.005, which corresponds to an energy 
scale of n ~ 10. 


KqT 1 eT)Xq cosh ? 7 Eq 

ef)Xq 2 E cosh ^ cosh ' 

(72) 

where r/ = and X q denotes an eigenenergy of the sys- 
tern. The intermediate steps of Eq. 72 are valid for any 
number of bands, but we have specialized to X q = ±E q 
for the particle-hole symmetric 2 -band case at the last 
step. Tr G q (r) is manifestly even in 77 , which guarantees 
that Tr G q (0) = Tr G q (/3). It can also be obtained from 
Eq. B 6 through direct simplification. 

Due to the presence of the energy scale set by T, there 
are two distinct regimes. In the (high energy) UV limit 
2 n < 2tt/3 the kinetic energy dominates the temperature, 
and we have cosh ~ e~ rEq + e^ r ~^ Eq : yielding the 

cosh —2“^ 

nice relation 

Tv C n (r) « Tr C n (r )| r=0 + Tr C n (0 - r)| T=0 . (73) 

This equation tells us that the nonzero temperature cor¬ 
relator above the energy scale of T is just the superpo¬ 
sition of two copies of the zero temperature correlator 
reflected about r = This is consistent with how cor¬ 
relation functions in BTZ geometry are obtained by a 
periodic quotient of those in AdS space[57]. Eq. 73 is 
completely general, since we haven’t used any particular 
form for the Hamiltonian. 

Inserting the result of Tr C n (r)\T=o from Eq. 49, we 
find that 

-log Tr C n {r) 

= - log (r -3 + (/3 - t ) —3 ) - 3(72 - 1) log 2 + log( 7 r) 

( 74 ) 


A, 


We can find the length scales R and £ r by comparing 
the maximal value of — logTr C n (r) with |^, r = We 
have 


'b/2 

Cx 


R v.-i 2 P 2 
Tr C ° Sh 
2R 2 p 
•7- lo g -r 
CT ^ 

2 R , 2 pf3 

V s TS 
2 R, 13 
{7 


(76) 


where we have used p = p n ~ approximate 

AdS geometry far from the horizon. Comparing this with 
— logTr C n {PI 2) = 3log ^ + log |, we obtain 


R 

Tr 


3 

2 


(77) 


which is exactly the same as in the T = 0 case. This must 
be true because the bulk geometry is still asymptotically 
AdS above the energy of T. However, the value of R is 
somewhat different due to the different functional forms 
of Eqs. 74 and 75, and takes the value 


R = {/T « 0.274 (78) 

V 7T 4 

In the IR limit below the energy scale of T, i.e. 2 n > 2i r/3, 
Eqs. 73 and hence 74 no longer hold because we must 
use the small f3E q approximation in the denominator of 
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Tr G q (r) = cosh m At r = f or p = 0, we have, to 2 nd 

cosh ^Z-SL z 

order in /3E q <C 1, 

Tr (r = 0/2) «^|l-^ j « 1 - y<? 2 (79) 
Hence, 


— logTr C(0/2) 

= — log Tr J dqG q (P/2)\W n (e iq )\ 2 
= - log (^j dq\W n (e iq )\ 2 -Lj dqq 2 \W n (e iq )\ 2 ^j 

~ - log (l-y</o|W n (e i9o )| 2 ) 


where we have used the fact that W n (e iq ) is sharply 
peaked at go = with magnitude Comparing 

Eq. 80 to the BTZ geodesic distance in Eq. 75: 


— logTr C(0/2) ~ El cosh 1 ^ « ^2 - l) (81) 

we arrive at 

pn0C K 1 + 9^^) (82) 

Indeed, p —>> b exponentially from above. This justifies 
the previous assumption that ^- 1<1 which is, among 
other implications, consistent with the fact that no non- 
negligible logarithmic subleading terms exist in Eq. 67. 
Physically, the agreement between the imaginary time 
distance of the bulk geometry and the BTZ black hole 
means that the rate of imaginary time correlator decay 
slows down exponentially as one goes deeper into the IR, 
which translates into the infinite redshift an outside ob¬ 
server sees for any physical process near the BTZ horizon. 


E. Critical boundary model with nonlinear 
dispersion vs. bulk Lifshitz black hole 

As a sequel to the previous subsection, we now consider 
a nonzero temperature boundary system with nonlinear 
dispersion in the long wavelength limit. As we have seen 
in the previous subsection, the energy-momentum disper¬ 
sion is sufficient for determining the decay properties of 
the correlators. Here we consider the simplest nonlinear 
critical dispersion 


E q ~ g 7 , for q « 0 (83) 

For higher g, E q should be regularized to a periodic func¬ 
tion. However, the details of the regularization do not af¬ 
fect the critical behavior as we have seen time and again, 


and Eq. 83 is sufficient as it stands. In this sense, Eq. 
83 subsumes the Dirac model (Eq. 12 ) with general M, 
which is merely an interpolation between a 7 = 1 and a 
7 = 2 Hamiltonian. 

The asymptotic behavior of the mutual information 
depends solely on the position of the singularities of 
G q (r) oc sech^|^-. They occur when f3E q = i(2l + 1)77 
l E Z, i.e. when 

E 2 = g 27 = — 7 t 2 T 2 ( 2 / + l ) 2 (84) 

In terms of z = e zq , they occur at zq = e *b 7rT ( 2 ^+ 1 )) 1/7 > 
Clearly, the l = 0 singularity has the largest magnitude 
within the unit circle, which is given by 


M = e' 


(TrTfi/^sin^ 


(85) 


From Eq. 27 and discussions surrounding Eqs. 65, we 
conclude that the mutual information between two points 
with angular separation of Ax sites decay like 

— log/ n (Ax) ~ 2 ( 7 rT) 1 / 7 ( 2 n Ax) sin — 

27 

= ( 7 rT) 1 / 7 sin (86) 

7T 27 

Similar to case (3) with linear dispersion, the AO depen¬ 
dence in the mutual information suggests that the cir¬ 
cumference approaches a finite value in the IR limit, so 
that there is an event horizon. However, the nonlinear 
dispersion leads to a different T dependence. Since the 
bulk geodesic distance oc — log/ n (Ax) oc T^A#, the 
black hole radius is b oc T 1 / 7 , different from the b oc T 
behavior of the BTZ black hole [58]. 

There are indeed black hole solutions to classical 
gravity with the property b oc T 1 / 7 , if one considers 
spacetimes with anisotropic scale invariance, i.e. with 
metric [59] 

r 2 7 c>2 2 

ds 2 =#- — —dt 2 + -o-dr 2 + — ^dx 2 (87) 

i? 2 ^ r 2 R 2 

invariant under the rescaling (£,x, r) —> (A 7 £, Ax, A _1 r), 
with 7 the dynamical critical exponent and R a length 
scale. 

If we include quadratic curvature tensor terms like Q 2 , 
Q , a /_or Ll a p /jLU Ll a/3 ^ l/ to the gravitational action[59, 
60], the resultant Einstein’s equations in the non 
Galilean-invariant spacetime will possess black hole so¬ 
lutions for certain ranges of parameters. Such solu¬ 
tions are known as Lifshitz Black Holes , which are well- 
studied [59-62] and proposed as possible gravity duals to 
Lifshitz fixed points in condensed matter physics. The 
explicit solutions of these black holes are known for cer¬ 
tain values of 7 , especially in the 2 + 1 dimensions rel¬ 
evant to our current context. For instance, the black 
hole metric for 7 = 3 and the gravitational action S = 
16 ^ G f d 3 x^/^g [f2 — 2T + 2 R 2 (fl a pQ a ^ — |$4 2 )], where 
f is the cosmological constant, is given by[59] 

ds2 = ~( x “ 7) df2 + + r2d ° 2 ( 88 ) 
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By examining the near-horizon geometry in Euclidean 
time, we explicitly find its Hawking temperature to be 
T = 2 ^ r4 , which agree with the horizon area we obtained 
from a boundary theory with cubic dispersion. 

The T oc b 1 dependence can also be expected from a 
simple counting argument. A system at a temperature T 
can be physically understood as one with states randomly 
distributed in a energy width of T. This randomness is 
quantified by the entanglement entropy S of the system 
with the thermal bath, most of which is carried by the 
IR region of the system. In the bulk system obtained 
through the EHM, the IR states carry maximal entropy 
per each site, so that the thermal entropy is proportional 
to the number of sites in the “stretched horizon”, which 
is proportional to the horizon area b. For a system with 
energy dispersion E oc g 7 , the momentum range that 
has energy below T is A k oc T 1 / 7 , so that the entropy 
S oc Ak oc T 1 / 7 . Consequently, b oc T 1 / 7 . 

The imaginary time bulk correlator behaves in a simi¬ 
lar way as that with linear dispersion (Case (3)). For the 
layers with energy scale above T, we still have, of course, 

Tr C n (r) = Tr [C n (r) + C n ((3 - r)] r=0 


X while the UV site now has a higher dimension y 3 , cor¬ 
responding to the I/E, JR, IR,UV and UV,UV sectors 
of the 1-dim case. More generally in D dimensions with 
a square lattice, one can map the 2 D sites in a cube to 
one IR site with Hilbert space dimension y and one UV 
site with dimension y 2 -1 . 

For free fermions systems, the mapping is equivalent 
to a wavelet transformation on the single-particle wave- 
functions, just like in the case with one spatial dimen¬ 
sion. The simplest higher-dimensional wavelet basis can 
be obtained via direct products of 1-dim wavelet bases. 
To define them, we first label the 1-dim wavelet functions 
as 




c^-^nUciz*- 1 ), v = i 

D^-'mUc^A v = 2 


(90) 


so that v = 1, 2 corresponding to the IR and UV wavelets 
in layer n, respectively. The //-dimensional wavelet func¬ 
tions can then be defined by 


D 


WZ"»- VD (z 1 , ...,z D ) = n W^( Zj ) (91) 

3 = 1 


which holds independently of the dispersion. In the IR 
limit with energy scale below T, the minimal value for 
Tr G q (r) at r = | still follows from Eq. 79 and 80, ex¬ 
cept that the energy eigenvalues are now A = q 1 . Hence 
we obtain a nontrivial (but still simple) nonlinear correc¬ 
tion 

(89) 

which is consistent with the metric in Eq. 88, which has 
the geodesic distance at r = | vanishing as a power of 
We still have p b exponentially as n increases, 
but at a different rate compared to the BTZ (Galilean- 
invariant) case. 


V. GENERALIZATION OF EHM TO HIGHER 
DIMENSIONS 

A. General setup 

When we generalize the boundary system to D spatial 
dimensions, the bulk system will contain D + 1 spatial 
dimensions, with a new emergent direction representing 
the energy scale. The boundary theory and bulk theory 
can be related by EHM in the same way as in the D = 1 
case. For example, with a boundary theory defined on 
the two-dimensional square lattice, a unitary mapping 
can be defined on four sites around a plaquette, which 
maps it to two sites representing the high energy and low 
energy degrees of freedom of the four sites. The mapping 
is illustrated in Fig. 5. If the Hilbert space dimension is 
X on each site, the output IR site should have dimension 


These 2 D wavefunctions for vj = 1,2 include one IR 
wavelet defined by Vj = 1, Vj and 2 D — 1 other wavelets 
that are regarded as UV degrees of freedom. The bulk 
correlators C^ u can be obtained from the boundary cor¬ 
relator G$(t) via this basis transform: 

C^(m,n 2 ,A(2-x),r) 

= E ^;(q)^ 2 (q)e iq ' A(2nx) G q (r) (92) 

q 

where we have denote the 2^ — 1 dimensional label of UV 
states viV 2 -..vd (with at least one Vj = 2) by p or v for 
simplicity. 



FIG. 5. Schematic picture of the EHM for two-dimensional 
boundary states. A unitary mapping (pink oval) is defined for 
four sites, each with a x-dimensional Hilbert space. It maps 
the DOFs from the four sites to two output sites, the IR site 
(red arrow) with dimension x and the UV site (blue arrow) 
with dimension y 3 . 

When there is K number of orbitals at each site, the 
correlation matrix is ( 2 D — 1)H x ( 2 D — 1)R For ana- 
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lyzing the asymptotic bulk geometry, it suffices to con¬ 
sider only the slowest decaying elements of which 

is determined by the lowest power of q in W%( q). Using 
the long wavelength asymptotic behavior C(e lQj ) « y/2 
and D{e iq i) ~ we see that among the 2 D — 1 UV 

wavelets W% lV2 " mVD ( q), the ones that play a leading role 
in the long wavelength correlation functions are those 
with only one Vj = 2 and all other Vk = 1, k ^ j. 
Therefore the asymptotic behavior of the bulk correlator 
is given by 


correlators u or v decay faster due to derivatives intro¬ 
duced by UV projectors D(e iq ). Their slowest-decaying 
elements involve D = 2 derivatives, i.e. d\, d\ or &XY- 
Hence 


M ~ M 


i 

(.x 2 + y 2 ) 3 / 2 


(95) 


I n ( Ax) ~ 6 log |Ax| + const. (96) 


C^(n, n, AX, r) ~ 2 D(n ~ 2) d Xj d Xk ]T e iq AX G q (r) 

q 

(93) 

where Xj = 2 n Xj, and j,k are the directions where a 
1 -dim UV wavelet function D is taken, i.e., Vj^ = 2. 

In the following, we shall examine the bulk geometries 
of various critical boundary systems at both zero and 
nonzero temperature, and highlight how they are differ¬ 
ent from those of (1 + l)-dim boundary systems. 


B. Critical boundary model at zero T 


Here, we shall examine in detail the decay properties 
of the bulk mutual information and bulk imaginary time 
correlator corresponding to a critical (D + l)-d boundary 
system. We will find that they describe a bulk geometry 
of a higher-dimensional AdS space, in close analogy to 
the (1 + l)-dim case described previously. Similar to one¬ 
dimensional case, we consider the Dirac model in (D+l)- 
dimensions [63-65]: 


* = £4 


D 


D 


sin ki 


M 


D=i 


D — cos ki 
i= l / 


Ck 

(94) 


where T°,T 2 are Hermitian Dirac matrices satisfying 
{T^jT 1 '} = for /q v = 0,1,...,D. For M close to 0, 
the lowest energy excitations of this system are centered 
around k = 0, where the EHM defined by wavelets in Eq. 
(90) and (91) correctly separates low energy and high en¬ 
ergy degrees of freedom. In the following, we study the 
behaviors of the correlation function and dual geometry 
along different directions. For simplicity, we shall only 
explicitly study the (2 + l)-dim case. 

Spatial (angular) directions. 

We build on the results of the decay of unequal spin 
propagators u, v for (1 + l)-dim, with two obvious ex¬ 
tensions mandated by Eq. 93: Firstly, we now need to 
perform a multi-dimensional sum over q and secondly, we 
need to decide which sequence of derivatives d\. produce 
the slowest decaying correlator. 

We recall that in the absence of a EHM transform, the 
critical correlator behaves like the inverse first power of 
distance, i.e. ~ , 1 Under the EHM transform, the 


which is an almost trivial generalization of the result in 
(l + l)-dim (Eq. 41). The undetermined constant defines 
the AdS radius of the corresponding AdS geometry, and 
is a complicated function of the full correlator involving 
u ^ v. Here, we have not been careful in keeping track of 
the powers of 2 n , and readers interested in doing so are 
invited to generalize the more rigorous derivation in Ap¬ 
pendix C. The apparent isotropy of Eq. 96 may not be 
exact due to the numerous approximations made. How¬ 
ever, any angular dependency should only manifest itself 
as a form factor in the correlators, with the leading log 
term in the mutual information remaining unaffected. 
Imaginary time direction. 

A critical D + 1-d boundary system also has power-law 
decaying imaginary time correlators, consistent with the 
interpretation of the bulk as a higher-dimensional AdS 
spacetime. Explicitly, 


/ on — 1 \ -^+2 

Tr C n (r) ~ — (97) 

V vft J 

with dimension-dependent critical exponent of D + 2. 
This is unlike that of the spatial correlators, which do 
not depend on D. Interpreted as a bulk geodesic dis¬ 
tance d r , we have 

= — logTr C n (r) ~ (D+ 2) [log(upT) — nlog2]+const. 

ST 

(98) 

which is proportional to the dimensionality. Physically, 
we can understand the origin of the D + 2 exponent as 
follows. Each spatial direction provides an additional 
dimension for the decay, and contributes a power of 
—— oc L . There has to be at least one direction where 

vpr vppr 

only the UV half of the degrees of freedom are selected, 
since a separation of energy scales is necessary for the 
EHM network. This direction contributes an additional 
power of 2 due to the gradient-like property of the UV 
projector D(z). The above statements are justified with 
more mathematical rigor in Appendix G 1, where the sub¬ 
leading powers in the decay are also explicitly evaluated. 


C. Boundary model with generic dispersion at 
nonzero temperature 

When an energy scale is introduced by the temperature 
T, the decay of correlators is dominated by the energy 
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scale which is independent of the EHM basis. Hence 
we can calculate the bulk correlators in a way similar 
to that of the (1 + l)-dim case, taking note only of the 
multidimensionality of q. 

Instead of just analyzing the Dirac model, we make our 
discussion more general by allowing our energy dispersion 
to take the following generic asymptotic form: 




n 


D 

E 




(99) 


i.e. with Vj = v and 7 j = 1 for j = 1 , 2 ,..., D. Let Ax be 
the displacement between two distant points within the 
same layer in the bulk. The mutual information decays 
like I xy ~ 8|^| 2 ~ 8|r?| 2 where, as shown in Appendix 

Dl, 

u ~ v ~ J e i2nq A5? tanh q 

( 102 ) 

Hence 


This form encompasses various physical scenarios, and 
reduces to the linear Dirac model in the simplest case of 
7 j = 1. When D = 1 and 7 > 1, Eq. 99 represents the 
nonlinear dispersions discussed previously. More inter¬ 
estingly, it can also describe semi-Dirac points character¬ 
ized by anisotropic dispersions, i.e. with D = 2, 71 = 1 
and 72 = 2 . Such dispersions have been observed in re¬ 
alistic systems involving ultrathin ( 001 ) VO 2 layers em¬ 
bedded in T 1 O 2 , which exhibit unusual electromagnetic 
properties [ 66 - 68 ]. 

According to Eq. 27 and subsection IV D, the correla¬ 
tors and hence mutual information decay exponentially 

o vp 

according to the complex root of cosh = 0 closest 
to the real axis. The roots are given by the values of q 
satisfying 


-log/ n (Af) 


2 n+1 7rT 

V 


I Ax I = 


LT 

~\ 


D 

E(A*i ) 2 (103) 

j 


This is a direct generalization of Eq. 66 for the (1+1)-dim 
critical Dirac model at temperature T, whose bulk geom¬ 
etry corresponds to that of a BTZ black hole horizon in 
the IR limit. Here, we have exactly the same asymptotic 
behavior, with the horizon having the same topology as 
the boundary system. When the latter is defined on a 
D-dimensional lattice with periodic boundary condition 
along all directions, the horizon is a D-dimensional torus 

j-iD 


2. ( 2 + 1)-dimensional anisotropic dispersion 


El = —7r 2 T 2 (2l + l) 2 , 

with l G Z, with the temperature T functioning as an 
imaginary gap. To find the exponential decay rate in 
direction j, we have to find \Im(qj)\, the imaginary part 
of the complex root q = qj of 


We now consider a generic anisotropic dispersion in 
D = 2 , and show that the event horizon can also become 
anisotropic. The dispersion is given by 

Eq = y+ 1 71 +1’!<72 72 , (104) 


v 2 q 2 ^ = -(to 2 + tt 2 T 2 (21 + l) 2 ) (100) 

where to 2 = X+ v 1 < h 1i denotes an effective mass from 
the momentum contributions from all the other direc¬ 
tions. Since T 2 and m 2 are both positive, we clearly 
choose l = 0 for qj to have the smallest imaginary part, 
i.e. slowest decay. While the decay rate also depends 
on m, the combination of momentum components giv¬ 
ing m = 0 yields the slowest decay rate \Im(qj)\. We 
can take m = 0 to be the dominant contribution to the 
overall decay rate hj, and the m > 0 contributions as 
the subleading corrections. This will be discussed explic¬ 
itly for the two cases below, with calculational details 
relegated to Appendix D 1 . 


1. Finite temperature Dirac fermions 


We first discuss the massless Dirac case with 


D 



( 101 ) 
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where j' = 1,2 for j = 2,1. It is mathematically tricky to 
obtain the asymptotic behavior of I xy for arbitrary Ax, 
when all components of x are not small. For our current 
purpose, it suffices to expand the asymptotic behavior 
about the limiting directions Ax = xe x and ye y . After a 
Gaussian integral computation detailed in Appendix D 1 , 
the mutual information at Ax = | A+|(cos + sin^e^) 
for <f> near 0 is approximately given by 


- log/n (Ax) 1^0 


2 n+1 


( • 

\ 271 



1/71 

X + 


71 ( 7 rT ) 2 - 1 ^ 1 + 71 yf\ 
2sin + v 2 x j 


= 2 n+1 \Ax\ sin — 
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1) 


4> 2 + o ftp 


(105) 


where aj = 

result holds near 
72 , a 2 and v 2 . 


(ttT) 




I /77 


——. An exactly analogous 

-J ! 2 7j 

(j) 1 , with 7 i,aq and replaced by 
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Notably, in the isotropic linear Dirac case where 7 j = 1 
and v\ = ^ 2 , OLj = 1 , the mutual information is manifestly 
asymptotically isotropic to third order by Eq. 105. This 
is despite the fact that the wavelet basis was constructed 
via tensor products of those of each direction and hence 
only possess four-fold rotation symmetry. 



FIG. 6. Plots of I xy ((j))/I xy (4> = 0) for a = 1 (almost hori¬ 
zontal line) and a = 2 (upper curved line) according to Eq. 
105, with higher order terms in qb 2 kept. We see that in the 
isotropic linear Dirac case with a = 1 , I xy oc |Ax| to a high 
degree of accuracy even away from ^ < 1. 


By contrast, when the dispersion acquires some non¬ 
linearity, 7 j > 1 and OLj 7 ^ 1 for some j and we expect 
the mutual information to be significantly anisotropic in 
a temperature dependent way. This is illustrated in Fig. 
6 , where the angular dependence of the mutual infor¬ 
mation is compared for a = 1 and 2. In any case, the 
factor 2 n+1 A|x| in Eq. (105) suggests that there is still 
a finite area (anisotropic) horizon in the IR limit, since 
the circumference of a closed circle around any periodic 


direction approaches a finite value while n —>> 00 . 


VI. CONCLUSION 

In this work, we have analytically studied the emergent 
bulk geometries of several different boundary systems 
through the EHM approach. In general, critical bound¬ 
ary systems at zero temperature correspond to scale in¬ 
variant bulk geometries. A spatial boundary appears in 
the infrared region when a mass scale is introduced by 
a nonzero mass. At nonzero temperature, a horizon ap¬ 
pears in the infrared region, which is distinguished from 
the spatial boundary by the infinite red-shift that can be 
observed in the behavior of correlation functions along 
the imaginary time direction. For critical boundary the¬ 
ories with different dynamical exponents in time spatial 
and temporal directions, the spatial geometry is similar 
but the space-time geometry depends on the dynamical 
exponent, at both zero temperature and nonzero temper¬ 
ature. Most of the above results qualitatively still hold 
true when the EHM is generalized to higher dimensions. 

A major open question concerns the dual geometry for 
a Fermi gas with finite charge density. The existence 
of a nonzero Fermi momentum makes it inappropriate 
to use the same EHM mapping defined here, since the 
long wavelength limit will no longer correspond to the 
low energy limit. A modified tensor network is required 
in order to describe the correct infrared physics near the 
Fermi surface. 
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Appendix A: Properties of the Mutual Information 

1. Relation to the single-particle correlators 

The mutual information as given by Eq. 5 contains 
the two-site entanglement entropy, which depends on the 
two-site correlator C xy via Eq. 9. 

It is a matrix of single-particle propagators [C xy \ij = 
(Pifij), where i, j G {x,y} and /3 x ,/3 y are the bulk anni¬ 


hilation operators. Below, we shall write 

C Cy V ) =C0 + V (Al) 

where C x and C y the on-site single-particle correlators, 
and C x - y is the (non-hermitian) single-particle propa¬ 
gator from site y to x. The off-diagonal contribution 
C x - y decays rapidly for large \x — y\, and I xy can be 
accurately approximated in that limit. Below, we write 
C xy = Co + V, where Co = C x 0 C y and V is a pertur¬ 
bation containing the off-diagonal parts C x - y and C\_ y : 


Ixy — S x 


Sy S xy 


= S X + S y + Tr ((Co + V) log(Co + V) + (1 - C 0 - V) log(l - C 0 - V)) 


^S x + S y + Ti ((Co + V)(log Co + Cq L V - -(C 0 -V) 2 ) 

+(1 - Co - F)(log(l - Co) + (1 - Co)- 1 ^ - |((1 - C 0 ) _1 y) 2 )) 

« (S x + S y - Tr (S x ® Sy )) + 2Tr V + Tr [V(log C 0 + log(l - C 0 ))] + ^Tr VCpV + ^Tr V(l - C 0 )~ l V 
= (S x + S V -S X - Sy) + ^Tr VCpV + -Tr V(1 - Cq^V 
» [V(C 0 (1 - Co))-^} 




Cx-y * i /tt ,, . Cy- X + (x y) 


'Cy(I-Cy) 


Tr [rt x Cy- x ] 


(A2) 


where C XJ C y are the single-particle onsite correlators, 
and C x - y is the single-particle propagator between the 
two different sites x and y. Note that we have implicitly 
assumed that Co and V commute while going from lines 
2 to 3, which holds in the IR limit. 


2. The mutual information in terms of the reduced 
density matrix 


The mutual information can also be understood as 
the Kullback-Liebler divergence [69] of the distributions 
described by p xy and p x p y , where p x is the single-site 
reduced-density matrix (RDM) and p xy is the two-site 
RDM. Writing the trace over one site as Tr x and two 


sites as Tr xy , we have 


Ccy — S x + Sy *S*xy 

= -Tr y p x log p x - Tr x p y log p y + Tr xy p xy log p xy 
= -Tr xy p xy log p x - Tr xy p y log p y + Tr xy p xy log p xy 

= Tr xy p xy log 

P*Py 


= E 


log^ 

. PxPyJ 


(A3) 


To maximize the mutual information / xy , we need p xy 
and p x py to be as different as possible. Since the log 
function drops steeply below unity, we especially want 
to avoid situations where the (square root of the) single¬ 
site RDM eigenvalue is large compared to that of the 
two-site RDM. Since p x = Tr y p xy , the above-mentioned 
situation is more likely when the RHS contains a large 
number of contributions from states belonging to differ¬ 
ent y. Hence we conclude that a maximal / xy must have 
minimal spread of p xy , i.e. have 2-particle states that 
are maximally entangled in the precise sense of Eq. A3. 
Note that this scenario represents a hypothetical opti¬ 
mum, and may not be realized the physical systems that 
we have discussed. 
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Appendix B: Correlators for two-band fermionic 
systems with particle-hole symmetry 


Here we derive the detailed results for two-band, 
particle-hole symmetric models studied in this work. 
They hold for generic two-band models at arbitrary 
chemical potential //, although we will ultimately use 
them only for the Dirac model at /i = 0. 

Due to particle-hole symmetry, the onsite correlator 
(projector) C x takes the following form in spin space: 


c x = (A A) 


a iA 
—iA a 


(Bl) 


where a and A are real, up to inconsequential corrections 
at nonzero temperatures. At zero chemical potential 
we always have a = The single-particle propagator is 
slightly more complicated: 


C X .y = OMJ) 


b iu 
iv b 


(B2) 


When fi — 0, the equal-spin propagator b = 0 and the 
unequal-spin propagators u,v are real. At nonzero 
u,v are not necessarily real but b is still real. \u\ ^ |u| in 
general, since C x - y does not have to be hermitian in spin 
space alone. When the onsite propagators C x = C y , as in 
the case when they are related by translation symmetry 
in the same bulk layer, we obtain upon substituting Eqs. 


G q (r) = ^cosh(r^ q )I+ ^ sinh(r.E q )^ ((l + 


In the limit of zero temperature (3 —>• oc, the equal-time 
correlator G q = G q (r = 0) tends to 

Gq G q | /j,=o + - Eq) ^sgn(/x)I + 

= 6(n- Eq)I + e(Eq- (sga(fl)I 

q (B7) 


The physical interpretation of the above is clear: When 
> \/jl\, G q is exactly the same as in the /i = 0 case. For 
< |/i|, G q projects identically to either both bands or 
none depending on the sign of /i. 


Bl and B2 into Eq. A2 

Ixy = 

[a(l — a) — A 2 ][2b 2 + (|^| 2 + |u| 2 )] + 2bA(l — 2 a)Re[u — v\ 
“ (^2_ a 2)( A 2_ (1 _ a )2) ' 

(B3) 

where a, A, 6, u, v are defined in Eqs. B1,B2 and Al. At 
zero chemical potential, b = 0 and a = |, and Eq. B3 
simplifies further to 


L xy \fjb=0 


4(H 2 + |^| 2 ) 

1-4 A 2 


(B4) 


Since the onsite contribution A does not vary with the 
spatial displacement x — y, the decay properties of I xy 
at /i = 0 depends almost entirely on u and v. 

From Eq. Bl, it is almost trivial to write down the 
expression of the single-site Entanglement entropy (EE): 

S x = -Tr (C x log C x + (I - C x ) log(I - C x )) 

= -2 Y, I A ± lo S A ± + (1 - A ±) !og(l - A±)] 

A± 

(B5) 

where X± = a=b \A\ are the eigenvalues of C x . When /i = 
0, a = \ and the EE is maximal at S™ ax = log 4 at A = 0. 
Indeed, without unequal-spin correlation (A = 0), we 
have no information about what is happening to the other 
spin state. When A is small, the maximal entropy is 
corrected according to S x ^ S™ ax — 2A 2 = log 4 — 2A 2 . 
S x vanishes when a = A = which produces a pure 
eigenstate oc |—) + |+). 

We now present the explicit forms of the correlators, 
denoted collectively as G q . Given a Hamiltonian fo q , 
G q (r) = + e^(^q-M))- 1 (Eq. 17) is explicitly 


sinh f3/i 


cosh /3E q + cosh /? /i 


hq sinh /3E q 
E q cosh /3E q + cosh /?/i 


(B6) 


For zero /i but nonzero temperature, we have 


G ( 


qlAt=o — 2^ ' 


TL tanh 4^ 

2E n 2 


(B8) 


Of course, This further reduces to the usual projection 
operator given by \ (l — ^-) at zero-temperature. 
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Appendix C: Derivation of the bulk mutual 
information for (1 + l)-dim critical boundary systems 
at T = 0 


We start from the following expression for the mutual 
information I xy between sites x and y (Eq. B4): 


L xy |/x=0 


4(M 2 + H 2 ) 

1-4 A 2 


8\u\ 2 ~ 8|^| 2 


(Cl) 


where, as introduced in the main text and the previous 
appendix, u and v are the unequal-spin single particle 
propagators, and A the unequal-spin onsite propagator 
which tends to zero beyond moderately large n. Below, 
we shall derive their asymptotic behavior in detail. 


1. Angular direction 


We evaluate Eq. 35 for a large angular interval of Ax 
sites by deforming the contour around the branch cut 
from z = 0 to z = oo. (looking like a tight-lipped Pac- 
man): 

u ~ v ~ ^ [ W*(z~ 1 )W n (z)z 2 Ax (y/z — Ve 27ri z) — 

2 Jo ^ 

= [ W:(z~ 1 )W n (z)z 2nAx ^ r dz 
Jo V z 

(W*{z- 1 )W n (z)z 2n )z 2n( - Ax - 1) - 1/2 dz 

1 

Q(z)z x dz (C2) 



where W n (z) = -J=D(z 2Tl 1 )]\ n j= l C(z 23 ') and 
C(z), D(z) = 1^. We have decomposed the integrand 

into a term Q(z) = W*{z~ 1 )W n {z)z 2,ri that does not have 
negative powers of z, and with X = 2 n (Ax — 1) — 1/2 
still very large. We next integrate by parts to get the 
asymptotic behavior of u ~ v: 


u 



Q( i) Q ; (i) , 

X + l (X + 1)(X + 2) 

0-0 +__+ ... 

(X + 1)(X + 2)(X + 3) 

Q"( i) 


(C3) 


Here, we have stopped at the 2nd derivative of Q, be¬ 
cause it is the lowest nonzero derivative at z = 1. 
Q( 1) = Q'( 1) = 0 because they each must contain at 
least one factor of W n { 1) or W*( 1), both of which are 
zero due to the presence of UV projectors D^z 2 ™ ). We 


also truncate off higher derivative terms as they contain 
higher powers of As one may expect, u or v depends 
exclusively on the behavior of the EHM basis at z = 1 or 
q = — i\ogz = 0, the IR point where criticality occurs. 

Substituting the explicit form of Q{z) and differentiat¬ 
ing, we obtain 


U v ~ 


Q"( i) 

x 3 


o2(n—1) i 

7r 

23n—4 

7T (2 n Ax) 3 
1 1 
167r ( Ax ) 3 



(C4) 


That n drops out is not a coincidence, but a manifesta¬ 
tion of scale invariance. Note that the presence of the 
critical point, which is a property of the Hamiltonian, 
only ensures that the first line of Eq. C2 will not evaluate 
to zero; the power-law decay rate is entirely determined 
by the analytic properties of the chosen EHM basis at 
that IR point. In this case, the mutual information be¬ 
haves asymptotically as 


Ixy — -^n(^^) ^ 


1 1 

32^ (Ax) 6 


(C5) 


2. Radial direction 


Here, we evaluate Eq. 36 for the mutual information 
between sites that are separated radially. When one of 
the layer lies at the UV boundary of the bulk (layer 1), 
u and v between layers 1 and n can be easily evaluated 
viz. 


u,v = —i <£ —WZ(z~ l )Wn(z)z^ 

7 1*1=1 2 

.-+(1 --- 2 " ') 2 




N 


=1 -2 2ttv / 2”+ t 1 ' 


z t 2 


= l 


dz 

1*1=1 z 2± * 
dz 


r (l — Zz 2 ™- +z 2n ) 


— i 


(2 z 2 "- 1 -z 2n )dz 


1*1=1 * 2± » 71*1=1 z 2± i 


(C6) 


The first integrand diverges when 2 —>> 0, so the contour 
should be inverted about \z\ = 1 and closed at infinity. 
The second one diverges when 2 —> 00 for n > 2, so the 
contour should be closed like a Pac-man. 

Performing the resultant real integrals analogously to 
Eq. C2, we obtain 


27rV2 n_1 


(C7) 
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67rV2 n_1 


(C8) 


The single-particle bulk propagator is not hermitian in 
spin-space, with u ^ v. The above are exact results, 
not asymptotic ones. However, exact results like these 
usually do not exist for more generic wavelet mappings. 
The resultant mutual information is 

4 y = /(l,n)~4(M 2 + M 2 ) = ^^ (C9) 


Appendix D: Calculational details for the bulk 
mutual information at nonzero temperature for 
critical systems in arbitrary number of dimensions 


In this appendix, we shall fill in the mathematical gaps 
in the derivations of various nonzero T results for gapless 
systems. Since the mutual information I xy between sites 
x and y is (Eq. B4), 


4(M 2 + M 2 ) 

1-4 A 2 


~ 8|^| 2 ~ 8|?;| 2 


(Dl) 


where A trivially approaches zero beyond the first few 
layers n, we will just need to find the asymptotic behavior 
of the unequal-spin single particle propagators (two-point 
functions) u rsJ V . 


1. Decay of nonzero T correlators for generic 
critical dispersions 

We consider the energy dispersion Eq. 99 




n 


D 

E 


4 + 


(D2) 


In the j th direction, the decay rate of the correlators u , v 
is given by the the imaginary part of the root q = <+ of 

v h 2lfj = -o 2 + 7f2T2 ) ( D3 ) 

nearest to the real axis, where m 2 = Yi^j v iQi 7i repre¬ 
sents an effective mass from the momentum contributions 
from all the other directions, m 2 <C T 2 is always satisfied 
for a bulk layer sufficiently deep in the IR, i.e. of suffi¬ 
ciently large n. In terms of the variable z = e iq = e iqj or 
q = qj = — used previously, the decay rate is given 

by — log \zo\, where zo is the root of 

+ 2 (oLj - 1 )z 2 + 1 = 0 (D4) 


tion and taking the imaginary part, 


\Im(qj(m))\ = -++ ((?rT) 2 +m 2 ) 


,2\l/(27i) 


ho 


D 


u, (i + I ( D5 ) 

hj 


where Uj = sin=(^) Tj ' and Vj = (2'y j (nT) 2 ) 1 . 

However, this is not the physical decay rate as it still 
depends on the other momenta. To obtain the physical 
decay rate, we integrate over the latter (taking j — 1 
without loss of generality): 


U rsj v < 


J e i 2 nqA *tanh++ D q 


n f<i<h 

i=2 J 


e -2 n \Im(q 1 (m))\Ax 1 e i2 n Ef> 2 Ax 3<l3 


a -2 n U 1 Ax 1 


D 

n 

J=2 

D 

n 

J=2 


dqe 


— {U 1 V 1 Vj2 Tl Ax 1 )q e i2 n Axjq 


_ e -2 n U 1 Ax 1 TT e 4 U lVl v2 2 -A*! 


-2 n u-\ Ax-\ _ -— _ - Axj - 

Z U\£±X\ 4U 1 V 1 Ax 1 2 v 2 


(D6) 


This is just Eq.105 for the general anisotropic case. For 
the isotropic linear Dirac case where Uj = ^ and Vj = 
2 ^t) 2 f° r T Eq. D6 nicely simplifies to 


D 


u ~ v ~ exp —2 n U\Ax\ — 


(A x j ) 2 


AUiViAxx <rri v, 

J>2 J 


. n ttT 2 n 

= exp —2 —Axi -7p 4- 

1 V 1 


v 2(ttT) 2 j>2 


y (A xjl 

f -7)2 


-2 n MX. a Xi 1+ 


= e 


Ef> 2 (Ax,y 


= e -2”^|AJ| 


(D7) 


This shows that the correlators and hence mutual infor¬ 
mation decay isotropically in A+ for the isotropic lin¬ 
ear Dirac model, at least in the neighborhood of e^, 
j = 1,2, ...,D. The nonlinearity of Eq. D6 is further 
explored in the main text around Fig. 6. 


Appendix E: Calculational details for the imaginary 
time correlator 


within the unit circle and closest to its boundary, with 
ocj — 2e 2 /7 (( 7 rT) 2 + m 2 ) 1 ^ 7 . Solving the above equa- 

v 3 


Here we shall present the full derivations of the more 
involved results on the imaginary time correlator C n (r). 
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We shall derive the results for an arbitrary chemical po- and 
tential /q so as to illustrate the interesting continuous 
crossover of C n (r) as a nonzero chemical potential is in¬ 
troduced to a critical system. 

For large r and fi > 0, Eq. B6 and B7 simplify to 


Ge q >h(t) 


e -^r e -rE q 


2 


Hq 


GE q <^(r) 



(El) Hence the full correlator given by 


C n (r) 


a -flT 


P dqlWnie^l 2 

J — 7T 


e(E q 


v)e~ TE « ( ! -^) +0(li-E q )e^ 



(E2) 


(E3) 


where 0 is the Heaviside function. We only have to 
care about the extreme IR region of this integral, since 
e~ EqT = e~ VF \ q \ T decays rapidly for large r in a critical 
system. 

To proceed further, we only need to understand the 
IR behavior of W n (z) = E^ 2 ”" 1 ) n?=i C^' 1 ). The 
chemical potential sets an energy scale that divides two 
qualitatively different regimes. When the layer n is be¬ 
low above the energy scale of /q i.e. n < n* where 2 n = 
pp 1 , W n (e iq ) is effectively dominated by a sharp peak at 


Qo = §ft < However, for n > n*, W n (e iq ) is governed 
by its analytic behavior near the IR point z = 1. Perturb¬ 
ing away from the IR point with z. = e h°+ A< d ~ 1 — iAxq 


|Wn(e iA9 )| 2 


\W n (l - *A)| 2 
W*(l+iAq)W n {l-iAq) 


1 

27r 



2 

| jD , (l)|2 2 2(„- 1 )(A g )2 


32^ (E4) 

_I 


Noting that the even part of ^ is the correlator 
simplifies to (with E v = (i) 


C n {r) 


e“^ r 2 3n 
32 t t 



0(q - v)e~ TV ™ (l - |a 2 ) + 9{v - q)e TV ™ (l + |a 2 ) 


(E5) 


This integral can be exactly solved in terms of incomplete 
Gamma functions. However, we just want to extract the 
relevant asymptotic behavior set by the scale fir. Since 
the correlator captures the IR behavior, it should remain 
invariant even if the upper limit of i r is replaced by an 
arbitrary q cu toff in the first term on the RHS. The fol¬ 
lowing formulae come in handy: j p q 1 e qr dq ~ for 
(it < 1 and for (it > 1. Also, J qcutoff q^ e - qr dq ~ 

P^jrr for gr < 1 and e ^ ^ for fir 1. 


text, (it = 0 and we always have 


Cn (l~) | /j, = 0 ~ 


rsj 


e — 2 3n 1 

167r VpT 3 

1 / 2 n \ 3 

167r \vpTJ 




1/4=0 


(E6) 


For the case of chemical potential discussed in the main 


The extension of this result to nonlinear dispersions will 
be discussed in Appendix F. Eq. E6 may also be used for 
the short-time behavior when g / 0, as long as /xr < 1. 
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Appendix F: Extension to the critical (1 + l)-dim 
Dirac Model with M ^ 1 and general discussion of 
criticality 

In the main text, we have focused on the M = 1 case 
of the critical (gapless) Dirac model 


Hoiradk ) = vp [sin kai + M( 1 - cos k)a 2 ] (FI) 

where M controls the relative weight between the sin k 
and 1 — cos k terms. These two terms have respectively 
linear and quadratic dispersions for small |fc|, and here 
we study the effects of their interplay. 

A simple plot reveals that the dispersion 

= y/ sin 2 k + M 2 ( 1 — cos k ) 2 looks almost per¬ 
fectly quadratic for M > 5. However, the short linear 
region near k = 0 is still expected to dominate the 
physics at the IR layers of the bulk. To see that this 
is indeed true, we explicitly calculate the order of the 
dispersion which is given by the derivative d l °^f w , 
where w = log k: 

d\ogE w _ e w (—M 2 + (-1 + M 2 ) cos [e™]) cot [^] 
dw —1 — M 2 + (—1 + M 2 ) cos [e w ] 

(F2) 

For small negative values of w = log fc, d l ° c f u Ew = 2(1 — 

M 2 e~ 2w ) = 2 (l - ^) « 2. But -S- 1 for large 

negative w. The transition region occurs at k c ~ as 
shown in Fig. 7. Note that a simple Taylor expansion will 
not reveal a quadratic dispersion, because it is unable to 
concentrate on an exponentially small IR region. 



FIG. 7. d l °d n Ew for TI as 2000. We see that the dispersion 
of Ek is linear (^f = 1) for w < — log M —7.6 or layer 
n ~ 14, but quadratic at momentum above that, till k ~ 0(1) 
where the dispersion levels off. 


a. Imaginary time correlator 


The zero-temperature asymptotics of the imaginary 
time correlation function are only dependent on the ex¬ 
treme IR behavior of the Hamiltonian. As evident in 
Eq. E3, a very large r in the imaginary time correlator 
suppresses contributions from all but the lowest energy 
regime. As such, a finite M should not change the long 
time behavior of C n (r). 

This is however not true for an infinite M, which will 
produce a purely quadratic dispersion since k c « = 0 . 

Let us write E & = vq k 1 for small &, where 7 = 2 here. 
Then Eq. E5 becomes 


o3n+l r 7r 

Tr C n {r) « J dqq 2 e-™°<>' 

= 27JRL 

167T j(vot) ^ 


(F3) 


Hence a purely nonlinear dispersion of order 7 affects 
the long time correlator by changing the exponent in the 
power-law decay from 3 to 3 / 7 . In our current context 
with the Dirac model, there exists a moderate imaginary 
time regime where the correlator decays like ~ r~ 2 . The 
duration of this regime becomes longer and longer as M 
increases, till it finally becomes infinitely long at M = 00 . 


b. Spatial correlators and criticality 

In general, the power law decay of the spatial correla¬ 
tors depends on the EHM basis, and cannot be changed 
unless the critical point becomes degenerate. To be pre¬ 
cise, the power law decay depends on the existence of a 
branch cut in the complexified correlator. 

In our two-band case, the complexified correlator 
which is introduced in Section IV just before Table I, has 
nonzero elements given by 

/sin k + iM( 1 — cos k) li(z — z~ x ) — iM(2 — (z + z -1 )) 

y sin k — iM( 1 — cos k) y i(z — z~ x ) + iM(2 — (z + £ -1 )) 

(F4) 

and its reciprocal. They have square-root branch points 
at z = 1, z = H 2+1 and . To compute the cor¬ 

relator, we perform the contour integral around \z\ = 1 
around the branch cut from z = 1 to z = , like what 

was done in Eq. C2. The correlator matrix elements are 
thus proportional to 

ft _ f j 7 

2 / 2 W* n {z- l )W n (z)z 2nA *y/F{z)- (F5) 

J 1-M 2 Z 

1 + M 2 


The shape of the dispersion discussed above affects the 
bulk geometry profoundly at nonzero temperature. The 

1 

emergent black hole radius behaves like 6 oc Tq when 
Ej~ oc k 1 . Hence M sets the critical temperature which 
separates the b oc T and b oc y/T regimes. 


(z $ ~ M l ) 

where F(z) = J-j — IXm* V ? or reciprocal. Via the 

V v - 1-M 2 ) 

same steps leading to Eq. C3, we will eventually find a 
^ ^3 decay in the spatial correlator, as long as the z = 1 
branch point is present. 
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The branch point at 2 = 1 may disappear when it com¬ 
bines with another branch point. In our case, it happens 
when M = 00 . Then F(z) becomes trivially equal to 
unity, and the correlator is identically zero. 

There are other more interesting degenerate cases 
where we end up with a spatial correlator that decays 
exponentially , even though the system is gapless. This 
happens when F(z) has singularities within the unit cir¬ 
cle, while the gapless point on the unit circle is not a 
branch point. An example is given by 

H(k) = sin 2 kai + (1 — cos&)<T 2 (F 6 ) 


where (letting vp = 1 for simplicity) 

£ 2 = +(2«j + l)(n — 1) (Q2) 

3 27r 

We next perform the integral in Eq. G1 iteratively, start¬ 
ing from the integral over q x : 

TYC(r) « f A 2 x A 2 y k 2 y ^J ky , Kx (r)dk y (G3) 

J —TV 

where Jk y ,K x is an effective massive (l + l)-dim correlator. 
For large r > — and k x = 1, it can be approximated by 


whose singularities in 2 = e lk occur at \i — l+iy/l + 2i\ = 
0.346, 1/0.346 and 1, with 1 being a double root that can¬ 
cels off in 1^. Hence its correlator decays like ~ 0.346 2Tla \ 
Physically, the gapless point at 2 = 1 is not critical be¬ 
cause the two bands touch but do not intersect. The ex¬ 
ponential decay arises from the effective mass scale due 
to the curvature of the dispersion. 

For two-band models, gapless points of even order are 
always noncritical (degenerate). However, such points 
may be critical if there are more than two bands. In 
general, an 7V-band gapless point will be noncritical if 
the order of its dispersion is a multiple of N. 


Appendix G: Derivation of results for 
higher-dimensional critical systems at T = 0 


1. Decay of the imaginary time correlator for 
general D 

It is instructive to first perform the derivation for a 
(2 + l)-dim boundary system. From Eq. E3, we have 

TvC n (r) = r r dq x dq y \W n (e iqx ) \ 2 \W n (e iqy )| 2 e _ ' EqT 

J — 7T J — 7T 

(Gl) 

where W n is the (1 + l)-dim bulk projector and Tr is a 
trace over the band indices, not the v indices (suppressed 
for now) labeling the 2 D — 1 bulk sectors containing var¬ 
ious combinations of one-dimensional holographic basis 
vectors. For large r, it suffices to consider the contri¬ 
butions close to the IR point q = 0 where W n (e iqj ) is 

maximal and F q ~ q 2 + q 2 . As explained in the main 

text, W n (e iq i) either behaves like a constant or is linear 
in qj near qj = 0, depending on whether the IR or UV 
projector is chosen. In this appendix, we shall derive the 
forms of all the terms in the correlator TrC^ 1 ^ 2 , and not 
just the dominant terms. 

Let us write the v index in binary form (fti, ft 2 , • ••, ftp), 
where ftj = 0,1 depending on whether the leading factor 
of W%(z) corresponds to an IR or UV projector. From 
Eq. E4 and the definition of the projectors in the main 
text, we know that \W n (e iqj )\ 2 ~ Cj4j* s f° r Qj < §?> 


Jk, Kx =i(r) = r e-V^-q^dq 

J —7T 

0 -ytk 2 +q 2 r q2k x ^ 


POO 

2 / e 
JO 


POO 

= 2 / e—(e 2 - k 2 Y 

Jk 

p 00 

/ e~ eT (e 2 - k 2 ) 

Jk 


\/e 2 — k 2 


de 


= 2 

« 2 


2 V ..-1 <^~k 2 ) 


de 


\J e 2 — k 2 

jTe-'V - t 2 )”*" 1 (< a - y) * 


= 2 f 


{e 2 -k 2 Y x + y(e 2 -)fc 2 )^- 1 


de 


= e 


-kr (2 + kr) 2 


0 — kr 


Qk x = 1 ( T k) 


(G4) 


where Q K is a 2 ft-th degree polynomial with constant 
term (2ft)!. The approximation from line 1 to 2 is ex¬ 
tremely accurate for large r, while that from line 4 to 5 
is valid for for extremely small k. This is the regime that 
contributes most to TrC n (r), because Jk, Kx is suppressed 
by at least like e~ kr where r is large. For small fc, the 
integrand does not decay fast, and indeed it is the regime 
where u^> k that contributes most to the integral. The 
other case, J kKx=0 (r), resist all known approaches of an¬ 
alytical approximation. However, it is obvious that it 
behaves asymptotically like 


Jk, K .=o(r) ~ *— (G5) 

from the relation (. d 2 - k 2 )J k:Kx=0 (r) = J k ^ x=1 (r), i.e. 
with Q Kx=0 (rk) ~ const. 

For a critical system in D +1 dimensions, we just have 
to replace E q by VF^^2f = i Q^-> an d substitute that into 

Eq. Gl. Now let’s define Tj (r) to be the integrand of 
Tr C n (r) with the first j dimensions integrated over. Our 
goal is to find the asymptotic behavior of Tjj{t). We 
have 


Ti(r) 


-Pi' 



(G 6 ) 
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where Pi = y 2}=i (Zj • We perform the integral over last 

variable using the same approximations (valid for large 
r) as in Eq. G4: 


valid for mr < 1 due to the approximations in Eq. G4. 
Hence 

Eq. G8 becomes 


Ti(r) 


i —1 


II f dqjChT j 

„_1 J — 7 T 


i=i 


dqCfq 2Ki e- p ' T 

J —TV 


TrC n (r)| 


0<mr<l 


-mr 2 ( 2 «j+l)(n-l) 

wn - 




= 2 C 


2 g 


n r d ^hY 3 

•_ 1 j—7T 


3 = 1 
-1 


/•oo __ 

/ dqq 2Ki e~V p n+^ T 

Jo 


3 =1 


(^t ) 2 


higher orders of 1 /r 

(G 10 ) 


and 


n r <%++ 

A —A 3 — 7T 


J =1 


— 2 Ci r 7 


d 


dr ) \ r 2Ki+1 


Q Ki (rPi_i) : 


F-i(r) 


0 ~-Pi-lT 


/y- 2 «i+l 


= 2 


z-l 


n^M r ^ r 


J=2 


r 2 (E15=2 k j)+ 




= 2 <( 2 ki)! n C| 

3=1 


3 =2 


—r 


dr 


r 2(52} = i 


OC 


id 


i 


j=i 


n 

3 = 1 


' J r 2 (Ej - =1 «j)+* 
— 1 \ 2 /^i + l 


(G7) 


On line 4, the tilde in Q K (—tJ^) denote normal ordering 
of the r and operators, i.e. and its products 

will have all r moved to the left of all When going 
from the sixth to seventh line, we note that each operator 
when acting on expressions of the 
form l/r At+1 , without incurring additional factors of 1/r. 

Hence 


TrC n (r)| mr >i ~ e mr xweak dependence on powers of 1/r 

The nonzero mass correlator is exponentially suppressed 
by e -mr , though for small mr we still see a subleading 
power law in r, albeit with a different power from that 
of the massless case. 


2. Nonuniversal properties of (2 + l)-dim critical 
model 

Here we illustrate how a different choice of model in 
(2 + l)-dim can affect certain quantities but not others. 
We consider the model 

H(qx,q v ) = d(q) ■ a 

= sinq x ai + smq y (j 2 + (cosg^ — cos q y )os 

(GH) 

which is also gapless at q = 0. Its correlator in momen¬ 
tum space is 

G q = \ (l - d(q) ■ 0-) 


1 cosq x —cosq y 

E~ q 

sin q x -\-i sin q y 
E~a 


sm q x — t sm q y 

E q 

. cosq x -cosq y 
1 E a 


I (G12) 


TrC n (r) ~ 


2 (2«jM)(n-l) 


1 D 
lpT D Pj 


(G8) 


, after restoring vp- Evidently, the leading terms occur 
when Kj = 1 for just one j, and are zero for the others. 
Hence, we have 


where E q = ^/2(1 — cos q x cos q y ) —> q 2 + q 2 = \q\ near 

criticality. When considered as a (1 + l)-dim correlator 
depending on q x (q y ) alone, it behaves like a massive cor¬ 
relator with mass q y (q x )j as can be seen from its poles at 
±icosh _1 sec q y (and vice versa for q x ++ q y ). 

Since the single-site correlator C x is given by 


Tr C n {r) 


(-) 

\v f t J 


D+2 


(G9) 


C r 


/:/: 


dq x dq y \ W n (e* 9x ) | 2 1(e i9y ) | 2 G q (G13) 


This is the main result for the imaginary time correlator 
in the multidimensional critical case. 

For the multidimensional massive case with a fixed 
mass m, Pi in Eq. G7 is replaced by Pi = 

yj‘m 2 + Y^j=i Qj- Hence Ti(r) also contains a mass and 
the third last line of Eq. G7 becomes 


2 i 


ik« 


3 = 1 



e mr 

7- 2 (P2j=i 


where \W n (e iqj )\ 2 is even about q = 0, we see that the 
off-diagonal components, being odd in q x or q yj must dis¬ 
appear. This is different from the 1 + 1-dim Dirac model, 
where ^(g), which is not odd in g, plays an important 
role in the decay of the correlator. In the current (2 + 1)- 
dim case, it is ds that controls the decay of the correlator. 

Juch as importantly, note that the nonconstant part 
of the diagonal terms are odd under the interchange 
Qx +>■ Qy Hence G q oc due to the symmetry be- 
twee the wavelet bases W n (e' lQx ) and W n (e iqy ). With the 
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off-diagonal part A (defined previously) vanishing rigor¬ 
ously, the single-site entropy is always maintained at ex¬ 
actly S x = log 4, the same universal limiting value in the 
(1 + l)-dim case. Evidently, the small n (UV) behavior 
of the entropy depends nonuniversally on the details of 
the model. 


Appendix H: Relationship between real and 
imaginary time correlators 

Our discussion of the EHM will not be complete with¬ 
out a proper discussion of the real time correlator, which 
is arguably of more direct physical significance. However, 
its oscillatory nature makes it unsuitable as a definition 
of bulk distance. Here, we shall discuss its mathematical 
and physical significance with the imaginary time corre¬ 
lator. 


1. Critical case with linear dispersion 

When there is a linearly dispersive critical point E q = 
vpq, Galilean invariance is restored and there is symme¬ 
try between space and time. Restricting ourselves again 
to (1 + 1)-dimensions, the bulk correlator within a layer 
is explicitly given by 

C(n,n, Ax, At) = J2\W n (q)\ 2 e i2nqAx G q (-iAt) 
q 

= Y,\ W n^\ 2ei2nqAXeiEqAtG 1 

q 

= \W n {q)\ 2 e iq ^ Ax+VFA ^G q 
q 

(HI) 

which depends symmetrically on 2 n Ax and vpAt. From 
Eq. C4, we deduce that the magnitude of the real time 
correlator behaves like 

1 / 2 n \ 3 

I C(n, -iAt)\ = u\ Ax=2 - nvpAt ~ — ( ( H2 ) 

which is identical to that of the imaginary time 
correlator [70]. This conclusion is consistent with the re¬ 


sult obtained by a naive Wick rotation, since extra factors 
of iin a power-law do not affect the decay behavior. 

When the dispersion is nonlinear, Galilean invariance 
is lost and the correct result cannot be simply obtained 
via Wick rotation, even if the system is still critical. 

2. Non-critical cases 

When a mass scale is present, the energy E q is bounded 
below by a value m, i.e. E q = m+e q where e q > 0. Hence, 
the real time bulk correlator 

c(m,n 2 ,0, At) =J2 W n 1 (QWn 2 (q)e iEqAt G g 
q 

= e imAt W^WnMe^Gg 
q 

(H3) 

acquires an oscillatory phase with frequency m, a result 
again consistent with Wick-rotating the exponential de¬ 
cay e~ rnr behavior of the imaginary time correlator. 

Appendix I: The behavior of geodesic distance 

Some results of this subsection and the next can also 
be found in Ref. 30. Here, we reproduce them for com¬ 
pleteness. 

1. AdS space 

Due to its remarkable symmetry, Anti-de-Sitter space 
can be embedded in a flat Minkowski spacetime one di¬ 
mension higher. As such, it inherits the simple metric 
structure of the latter, which yields simple expressions 
for geodesic distances. 

For brevity, lets only consider the Euclidean (2+l)-dim 
AdS space, since its higher dimensional analogues will 
give rise to very similar expressions. We parametrize the 
space by coordinates w = (p, #, r), and embed it in a 4-d 
Minkowski spacetime as the locus of X a x br rj a i ) = X • X = 
R 2 , where R is the AdS radius, p = diag( 1, —1, —1, —1) 
and 


X fl (w)= vV + i? 2 cosh 


P2 I JA- 
n T 47r 2 


, \Jp 2 + R 2 sinh 


R2 I JA 
n + 4 tt 2 


= , pcosO, psmO 


(II) 


Here r appears with a rescaling factor of , 1 == , in- 

X+E 

stead of the more conventional i • This is to ensure that 


the rescaled AdS metric 

1 + fe dr 2 + -AE + fd# dr 2 + p 2 de 2 (I2) 

1 L 2 1 I P 2 r r v 7 

1 + 47 X 2 R 2 1 + R2 
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is 0(2)-invariant at p = -R R. From well-known 
properties of the Minkowski metric, the geodesic distance 
between points wi and is given by 


d \2 = R cosh 1 




(13) 


From Eq. 13, we find the geodesic distance corresponding 
to an angular displacement to be 


<C” = flcoah-' 


= Bcosh -’(l + % ' ,M ' 


sin 


= 2iilog^ 


(14) 


is given by 


d T = R cosh 1 
~ R cosh -1 


((g + 1 ) c " h tf+ 


R cosh 1 


(P 2 (27r) 2 T 2 \ 

Vi? 2 2(L 2 +47r 2 i? 2 )y 


~ 2Rlog 


2irpr 

RL 


R 2 


( 16 ) 


where we have used p R while going from line 1 to 2, 
and t L from line 2 to 3 and pr RL from line 3 to 
4. 


2. Geodesics near a black hole horizon 

We start from a metric 


for p^> R. In the scale-invariant case where Ax = pAO , 
and we also have d~ 2 R log . 

For a radial displacement with A p = |p 2 — Pi|, we have 


d™ n = R cosh 


-l 


' yfW+fflW + pfi-pin 


-1 ( Pi P2 


i?cosh 


R R 

4 (ei 

,2 VP2 


R 2 

R 2 


R cosh- 1 ' ' ; 1 • ", II ", -1 


2 Pi 


El 

Pi 


1 P 1 
log — 

P2 


= R 
= ^|A(logp)|, 


R 2 


2 Pi 


(15) 


also for pi , p 2 ^ R . 

The geodesic distance in the imaginary time direction 


ds 2 = V(r)dT 2 + —— + r 2 dfl 2 (17) 

V(r) 

with V(ro) = 0, which admits a horizon at r — vq. 
This horizon is a null surface within which ds 2 = 0. 
Intuitively, the geodesics between two points infinitesi¬ 
mally close to the horizon must necessarily wrap around 
the horizon (i.e. with with r being constant), since the 
cost of radial displacements A r diverges to infinity at 
the horizon. Hence we always have As ~ r^AO near a 
horizon. 

To explore the near-horizon geometry more rigorously, 
we switch to Rindler coordinates valid near the horizon 
r = ro. We define 


p = 2sgn(r - r 0 ) 


\r - ro | 

innOI 


so that, to order 0(p 2 ), the metric becomes 


(18) 


ds 2 \ R =p 2 




dT 2 +11 — 


V 




dp + r 0 ± 


|V(r 0 )| p 2 


dn 2 


(19) 


where T = — ^— r and the ± sign refers to the region 
immediately outside and inside the horizon respectively. 

Upon dropping the 2nd order terms, we recover the 
usual Rindler metric which just describes a plane with 
T taking the role of the polar angle. To avoid a conical 
singularity, we require that T has period 27r, i.e. that the 
period of r and thus /3 = T -1 is of the value y 

We now specialize to an example most relevant to the 
main text, which is the near-horizon geometry of a (2+1)- 
dim BTZ black hole. The horizon occurs at r = 6, with 
another parameter R setting the overall length scale. The 


horizon must also occur at a D + 1-th order zero of V ( r ). 
Hence we have 


3 _ r3 

W (110) 


with V'(b) = and V"(b) = 0. That the second deriva¬ 
tive is identically zero is unique to D = 2 dimensions. 
After some tedious derivation, the geodesic distance at 
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Rindler radius po between two points separated by A 6 is 
AA 

R 


where b is the horizon radius and 




: cosh 


r/ z 


y/l + T] 2 cosh 



bA6 


V 


3 po , 3 bA6 /T _ N 

2^Wo^ ( I12 ) 


2 R 


/ J A moment of calculation reveals that the above reduces to 
(I11)AA = bAO as we approach the horizon where rj oc po —;> 0. 



